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The  turbulence  phenomena  associated  with  the  mixing  of  fuel  and  air  in  a  supersonic 
combustion  ramjet  are  studied.  The  turbulence  generated  by  coaxial  and  normal  jet  fuel 
injection  is  studied  by  analyzing  a  control  volume  in  which  fuel  and  air  are  mixed.  These  two 
modes  of  fuel  injection  are  eximined  for  various  initial  conditions  in  order  to  compare  their 
respective  characteristics  and  determine  the  design  advantages  of  each  mode. 

The  results  of  this  analysis  show  that  large  turbulence  energies  can  be  generated,  much 
larger  than  the  energy  added  t<y  the  jet.  Thus  one  draws  a  mental  picture  of  the  fuel  jet 
exciting  the  flow  into  a  turbulent  state.  This  turbulence  energy  if  frozen  into  the  flow  causes 
losses  In  thxuat  and  specific  impulse.  These  losses  are  calculated  as  flight  speed  varies  for 
both  modes  of  fuel  injection.  A  universal  plot  for  a  diatomic  gas  Y  -  1.4  was  determined  for 
which  specific  impulse  losses  are  found  to  be  a  function  of  relative  turbulence  intensity  and 
flight  speeds  These  losses  are  large  but  not  catastrophic  as  previously  reported  by  Dr.  James 
Switheabank  of  the  University  of  Sheffield,  England. 

Experimental  results  from  a  compressible  turbulent  shear  flow  experiment,  the  decay  of 
a  supersonic  free  jet,  having  important  theoretical  implications  and  are  used  in  formulating  a 
new  theoretical  model  for  turbulent  shear  flow,  A  basic  relationship  between  the  local  "eddy 
viscosity"  and  the  square  root  of  the  local  turbulent  kinetic  energy  has  been  found  for  in¬ 
compressible  and  supersonic  jets.  This  result  plus  a  correlation  of  the  turbulent  energy 
dissipation  rate  allows  the  formulation  of  a  turbulent  flow  model  that  includes  calculation  of 
the  turbulent  kinetic  '  -  field  as  well  as  the  mean  velocity,  mean  temperature,  and  mean 
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13.  ABSTRACT  (C'mtd) 


A  fuel  injector  resigned  to  maximize  the  turbulent  kinetic  energy  was  tested  to 
verify  the  existing  n  ixing  rate  correlations.  This  fuel  injector,  a  normal  jet,  of  high 
Mach  number  and  re  :tangular  cross  section  exhibited  slightly  larger  penetration 
distances  than  predicted  by  a  recent  correlation  developed  by  Povinelli  at  the  NASA 
Lewis  Flight  Propulsion  Laboratory.  The  mixing  rate,  as  measured  by  decay  of  the 
maximum  eoneeutra:  ion.  was  nearly  an  order  of  magnitude  faster  that*  predicted  by  the 
correlation  due  to  H*jnry  at  the  NA*IA  Langley  Research  Center.  Further*,  it  is 
postulated  that  the  mixing  length  A  this  injector  may  be  independent  of  initial  conditions 
The  results  of  this  Si.udy  are  directly  applicable  to  the  design  of  supersonic  combustion 
ramjet  engines. 
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ABSTRACT 


The  turbulence  phenomena  associated  with  the  mixing  of  fuel  and  air  in  a 
supersonic  combustion  ramjet  are  studied.  The  turbulence  generated  by  coaxial 
and  normal  jet  fuel  injection  is  studied  by  analyzing  a  control  volume  in  which 
fuel  and  air  are  mixed.  These  two  modes  of  fuel  injection  are  examined  for 
various  initial  conditions  in  order  to  compare  their  respective  characteristics 
and  determine  the  design  advantages  of  each  mode. 


The  results  of  this  analysis  show  that  large  turbulence  energies  can  be 
generated,  much  larger  than  the  energy  added  by  the  jet.  Thus  one  draws  a 
mental  picture  of  the  fuel  jet.  exciting  the  flow  into  a  turbulent  state.  This 
turbulence  energy  if  frozen  into  the  flow  causes  losses  in  thrust  and  specific 
impulse.  These  losses  are  calculated  i  s  flight  speed  varies  for  both  medes  of 
fuel  injection.  A  universal  plot  tor  a  dia  omic  gas  y-  1.*.  was  determined  for 
which  specific  impulse  losses  are  found  to  be  a  function  of  relative  turbulence 
intensity  and  flight  speed.  These  losses  are  large  but  not  catastropic  as 
previously  reported  by  Dr.  James  Swithenbank  of  the  University  of  Sheffield, 
England. 

Experimental  results  from  a  compressible  turbulent  shear  flow  experiment, 
the  decay  of  a  supersonic  free  jet,  have  important  theoretical  implications  and 
are  used  in  formulating  a  new  theoretical  model  for  turbulent  shear  flew.  A 
basic  relationship  between  the  local  "  eddy  viscosity"  and  the  square  root  of  the 
local  turbulent  kinetic  energy  has  been  found  for  incompressible  and  supersonic 
jets.  This  result  plus  a  correlation  oi  the  turbulent  energy  dissipation  rate 
allows  the  formulation  of  a  turbulent  flow  model  that  includes  calculation  of  the 
turbulent  kinetic  energy  field  as  well  as  the  mean  velocity,  mean  temperature, 
and  mean  concentration  fields. 


A  fuel  injector  designed  to  maximize  the  turbulent  kinet.j  energy  was  tested 
to  verify  the  existing  mixing  rate  correlations.  This  fuel  injector,  a  normal 
jet,  of  high  Mach  number  and  rectangular  cross  section  exhibited  slightly  larger 
penetration  distances  than  predicted  by  a  recent  correlation  developed  by 
Povinelii  at  the  NASA  Lewis  Flight  Propulsion  Laboratory.  The  mixing  rate, 
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as  measured  by  decay  of  the  maximum  concentration,  was  nearly  an  order  of 
magnitude  faster  than  predicted  by  the  correlation  due  to  Henry  at  the  NASA 
Langley  Research  Center.  Further.,  it  is  postulated  that  the  mixing  length  of 
this  injector  may  be  independent  of  initial  conditions.  The  results  of  this  study 
are  directly  applicable  to  the  design  of  supersonic  combustion  ramjet  engines. 
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SYMBOLS  (Contd) 


normal  coordinate  mesh  indice  for  numerical  analysis  of  differential 
equations 


mass  flow  rate 


downstream  mesh  irdice  for  numerical  analysis  of  differential 
equations 


number  of  eddies  crossing  station  2  per  unit  time  through  area  L 
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number  of  eddies  containing  air 

number  of  eddies  containing  fuel 

outer  normal  unit  vector  to  surface  element 

outer  normal  to  control  volume  at  station  1 

outer  normal  to  control  volume  at  station  2 

3tatic  pressure 

turbulence  static  pressure 
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static  pressure  outride  mixing  zone 

total  pressure 

total  pressure  of  fuel  at  injection  station 
total  pressure  of  air  at  station  1 


staiifc  pressure  at  station  1 
e^rrace  static  pressure 
time  avera  ;  total  pressure 


time  average  total  pressure  in  turbulent  flow  with  frozen 
turbulence  energy 


mean  static  pressure  at  station  3 
dynamic  pressure  of  fuel  at  injection  plane 
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SYMBOLS  (Contd) 


dynamic  pressure  cf  air  at  station  1 
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turbulent  heat  flux  vector 


specific  gas  constant 

radial  coordinate  in  cylindrical  polar  coordinates 

surface  of  integration 

effective  Schmidt  number 

static  temperature 

time  period  for  time  averaging 

Eulerian  integral  time  scale  of  turbulence 

turbulence  temperature 

effective  total  temperature  in  frozen  turbulent  flow 


total  temperature 
integration  time  for  i-th  eddy 
Xj  coordinate  velocity  of  U 

r  coordinate  velocity  in  cylindrical  polar  coordinates 

velocity  of  i-th  eddy 

air  velocity  at  station  1 

average  velocity  at  station  2 

a vet age  vd  .-city  at-  Btstion  *= 

time  average  of  velocity 

fluctuating  velocity 

root  mean  square  turbulence  velocity 

fluctuating  velocity  of  i-th  eddy 
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SYMBOLS  (Contd) 


U 


V 

V 

V 

V 
v' 
w 
w 

Wmox 

W 

w' 

X 

*0 


x„x2,x3 


Yf 

Yi 

Yo 

Y, 


z 


vector  velocity  (U,  V,  W) 
volume  of  integration 
X2  coordinate  velocity 

d  coot  iinate  velocity  in  cylindrical  polar  coordinates 

I 

mass  average  velocity 
fluctuating  component  of  velocity  V 
X3  coordinate  velocity 

z  coordinate  velocity  in  cylindrical  polar  coordinates 

largest  value  of  W  at  a  given  station 

time  average  velocity  in  cylindrical  polar  coordinates 

fluctuating  component  of  velocity  W 

variable  of  polynomial 

last  distance  from  normal  injector  for  which  concentration  was 
100% 

rectangular  coordinates 

mass  fraction  of  fuel  in  mixture 

mass  fraction  of  species  i 

time  average  mass  fraction  of  species  i 

penetration  height  based  on  zero  concentration  point 

mass  fraction  of  air  in  mixture 

transformed  str  •unwise  coordinate 

streamwise  coordinate 

GREEK  SYMBOLS 
ratio  of  specific  heat 

momentum  parameter,  1  for  coaxial  mixing,  0  for  normal 
injection 
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GREEK  SYMBOLS  (Contd) 


turbulera  stress  tensor 

shear  part  of  turbulent  stress  tensor 

shear  stress  tensor  of  turbulent  velocity  fluctuation 


nondimensionaJ  enthalpy 
represents  any  variable 

transformed  radial  coordinate  in  Von  Mises  plane 
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PREFACE 

The  prediction  of  turbulent  mixing  rates  has  always  been  one  of  the  darkest 
secrets  in  the  black  art  of  propulsion.  Because  of  its  importance  and  lack  of 
understanding  of  its  mechanism,this  phenomenon  has  been  restudied.  A  new 
formulation  is  developed  in  this  dissertation  based  on  new  data  obtained  from  a 
simple  turbuient  shear  flow,  the  supersonic  free  jet.  This  formulation  grew 
from  the  need  to  explain  the  static  pressure  defect,  measured  on  the  centerline 
of  the  free  jet  extending  from  the  end  of  the  potential  core  downstream  to  the 
last  measurement  station.  This  defect  was  measured  for  all  initial  Mach 
numbers,  total  pressures,  and  total  temperatures  tested  in  the  laboratory. 

This  pressure  defect  is  so  large,  in  the  high  Mach  number  region  of  the  jet, 
that  it  must  be  included  in  the  dynamic  of  the  flow.  This  static  pressure  defect 
is  now  believed  to  be  a  direct  measurement  jf  the  radial  component  of  the 
Reynolds  stresses  at  that  point  in  the  flow. 

I  wish  to  express  my  appreciation  to  Professor  Rudolf  Edse  for  the  friend- 
ship,  guidance,  and  understanding  he  has  given  me  as  my  advisor  since  '.he  start 
of  my  academic  work  at  The  Ohio  State  University. 

This  work  was  accomplished  at  the  Air  Force  Aero  Propulsion  Laboratory 
as  part  of  Project  3012,  Special  Ramjets.  I  also  wish  to  express  appreciation 
to  the  following  men  of  the  Ramjet  Division  for  their  contributions  to  this  work: 
To  Mr.  John  T.  Hojnacki,  friend  and  confidant,  who  taught  me  much  about 
modern  electronic  instrumentation  systems  and  was  always  there  when  extra 
help  was  needed.  To  Mr.  Wayne  A.  Zwart,  formerly  of  the  .Ramjet  Division, 
who  did  a  competent  job  of  mechanical  design  of  the  water  cooled  free  jet 
combustor  and  nozzle. 


To  Mr.  Kenneth  G.  Schwartzkopf  who  programmed  the  transformation  of 
the  data  acquisition  system  tapes  into  dimensional  flow  variables;  a  program 
that  I  wanted  changed  constantly.  To  the  technicians  of  the  Laboratory,  Mr. 
Robert  A.  Schelenz,  Mr.  Dante  Cincerolli,  and  Mr.  Ashton  Sayre  who 
cooperated  fully  with  me  to  put  forth  our  best  efforts  on  the  success'::*  -xperi- 
ment s  and  also  the  many  frustrating  and  unsuccessful  experiments  not  reported. 
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I  also  want  Co  express  my  appreciation  to  Dr.  Harold  Wright  of  the  Air  Force- 
Institute  of  Technology  who  supported  some  preliminary  hot  wire  anemometer 
experiments  in  the  Mechanical  engineering  Lab. 

Finally,  ?  want  to  thank  Susan,  for  giving  me  five  beautiful  babies,  our  most 
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SECTION  I 
INTRODUCTION 
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X,  1  BACKGROUND 

The  development  of  the  supersonic  combustion  ramjet  engine  for  hypersonic 
flight  propulsion  has  necessitated  much  research  on  the  processes  of  fuel 
injection  and  mixing  because  of  their  great  influence  on  engine  performance  and 
design.  Incomplete  mixing  is  obviously  undesirable  since  unmixed  fuel  cannot 
be  burned.  Long  mixing  lengths  not  only  increase  engine  weight  and  cooling 
requirements  but  also  drag.  Therefore,  rapid  mixing  is  desirable.  The  fuel 
injection  mode  strongly  governs  total  pressure  losses  due  to  mixing  so  that  total 
engine  performance  may  be  reduced  greatly  by  mixing.  The  mixing  flow  field  in 
the  supersonic  combustion  ramjet  is  highly  turbulent  and  thus  mixing  rates  are 
determined  by  the  turbulent  transport  of  mass  and  momentum  rather  than 
molecular  transport.  Turbulent  transport  results  from  the  time  correlation  of 
the  unsteady  velocity  and  concentration  fluctuations  in  the  flow  field.  The 
turbulence  kinetic  energy  produced  by  mixing  can  be  a  large  fraction  of  the 
stream  total  energy.  In  recent  publications  (References  1  and  2),  Swithenbank 
has  calculated  the  turbulence  generated  by  coaxial  injection  of  fuel  in  the 
combustor.  His  analysis  shows  that  large  losses  may  be  incurred  by  the 
turbulence  energy  production  due  to  mixing.  He  also  shows  that  this  loss 
increases  with  flight  speed,  and  that  it  ultimately  becomes  the  primary  factor 
limiting  the  top  speed  of  the  supersonic  combination  ramjet.  The  validity  of 
Swithenbank* s  results  are  open  to  criticism  since  his  analysis  is  based  mainly 
on  intuitive  iormula  rather  than  rigorous  theoretical  development.  Both 
experimental  and  analytical  studies  of  gaseous  mixing  have  been  used  for  engine 
design.  The  most  extensive  work  has  been  directed  at  the  coaxial  jet. 

Analytical  models  loying  many  eddy  viscosity  models  have  been  developed. 

The  correct  formulation  of  an  eddy  viscosity  has  been  of  prime  importance 
in  many  researches.  The  eddy  viscosity  models  developed  in  this  country  are 
generally  minor  modifications  of  the  Prandtl  velocity  defect  formula.  These 
formulas  relate  the  eddy  viscosity  to  the  mean  flow  velocity,  however,  rather 
than  to  the.  turbulence  velocities.  This  appears  to  be  a  basic  limitation  to 
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these  formulations  and  the  number  of  models  developed  attests  to  the  lack  of 
confidence  in  their  general  applicability.  This  Jack  of  confidence  stems  from 
the  fact,  that  no  fundamental  theory  exists  to  explain  these  models;  however,  an 
attempt  to  formulate  a  uniform  model  has  been  proposed  by  Shetz  (Reference  3). 
Recently,  a  new  approach  lias  teen  proposed  by  Bradshaw  and  Ferriss 
(Reference  4)  which  does  away  with  the  eddy  viscosity  concept  and  relates  the 
turbulent  shear  stress  directly  to  the  turbulent  kinetic  energy.  This  method 
has  teen  extended  to  coaxial  mixing  at  low  subsonic  Mach  numbers  by  Lee  and 
Harsha  (Reference  5).  This  method  has  the  disadvantage  that  turbulent  shear 
stress  can  only  be  zero  when  the  turbulent  kinetic  energy  is  zero  unless  ad  hoe 
assumptions  are  used  to  avoid  this  difficulty,  in  addition  this  model  appears  to 
reverse  the  roles  of  shear  stress  and  energy  in  the  dynamics  of  the  generation 
of  turbulence. 

Normal  fuel  injection  has  also  been  investigated  extensively  and  mathematical 
models  for  penetration  have  teen  developed  (Reference  6);  however,  no 
analytical  solutions  of  the  mixing  problem  have  been  obtained  due  to  the 
mathematical  complexity  of  the  three  dimensional  flow  field.  An  important 
correlation  of  mixing  rate  for  normal  jets  at  various  initial  conditions  has  been 
obtained  by  Henry  (Reference  7).  Development  of  practical  mixing  systems 
depends  primarily  on  empirical  results  for  both  coaxial  and  normal  injection. 
These  two  modes  of  injection  involve  very  different  mixing  characteristics  and 
pressure  losses.  Figure  1  in  the  appendix  shows  the  structure  of  turbulence 
obtained  by  a  spark  Schlieren  photograph.  The  highly  turbulent  region  in  the 
center  of  the  photograph  is  the  turbulence  generated  by  the  shear  flow  of  a 
superson*c  free  jet.  It  is  evident  that  any  theoretical  calculation  of  turbulent 
flows  should  include  a  calculation  of  the  turbulent  energy  field  also;  this  will  be 
the  main  effort  of  this  dissertation. 

1.2  PURPOSE 

The  purpose  of  this  dissertation  is  to  investigate  the  influence  of  the  initial 
conditions  in  supersonic  combustors  on  total  pressure  loss,  the  production  of 
turbulence  and  performance  loss,  and  then  to  determine  the  basic  relationship 
between  turbulent  kinetic  energy  and  turbulent  transport  rales. 
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SECTION  II 

THEORETICAL  RELATIONSHIPS  FOR  TURBULENCE 
PRODUCTION  AND  PERFORMANCE  LOSSES 


2. 1  DISCUSSION 

The  purpose  ot  this  section  is  to  rigorously  formulate  the  mathematical 
analysis  in  order  to  first,  calculate  the  maximum  turbulent  kinetic  energy 
produced  in  a  given  mixing  flow,  and,  second,  to  calculate  the  effect  of  this 
turbulence  on  the  performance  of  the  supersonic  ramjet  engine. 

This  analysis  is  formulated  for  both  the  c  axial  and  normal  fuel  injection 
mode  so  that  the  importance  of  geometric  boundary  conditions  can  be  assessed. 
The  calculations  are  also  done  for  incompressible  and  compressible  flow  so  that 
the  effects;  of  compressibility  can  be  determined  as  well  as  the  effects  of  the 
dynamic  initial  conditions;  Mach  number,  molec-ular  weight,  total  temperature, 
and  toial  pressure.  The  total  pressure  losses  will  be  studied  to  determine  if 
turbuleuce  intensity  can  be  varied  independently  of  total  pressure  loss.  The 
theoretical  results  of  this  section  are  then  used  to  determine  the  single  element 
injector  design  with  the  highest  turbulence  intensity. 

2.2  ANALYSIS  OF  INTERNAL  FLOW 

The  production  entropy  can  occur  by  reversible  and  irreversible  processes. 
The  irrevetoible  processes  are  caused  by  gradients  of  concentration  arid 
velocity  driving  the  molecular  transport  of  mass,  energy,  and  momentum.  We 
will  assume  that  this  entropy  production  is  the  result  of  mixing  and  the 
dissipation  of  turbulent  energy.  Thus  we  can  calculate  the  total  turbulence 
intensity  produced  in  a  given  flow  by  considering  flows  without  entropy  change. 

The  one-dimensional  or  integral  equations  of  motion  for  turbulent  flow  can 
be  easily  developed  by  foliowing  the  procedure  of  defining  the  instantaneous  field 
variables  to  lie  the  sum  of  the  time  mean  and  fluctuating  part  and  integrating  the 
equations  of  motion.  This  integration  defines  the  appropriate  average  value  at 
any  streamwiae  location. 
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Since  thif?  integration  is  both  a  time  and  space  integration,  the  equations  are 
mean  conservation  equations  and  they  include  the  Reynolds  stresses. 
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Sketch  1.  Control  Volume  for  Constant  Area  Mixing 
The  continuity  equation  is 

r  a/>  ,  r  ^  „ 


/  — ^  dV  +/  RU  •'n  dS  =  0  ( 

v  dt  s 

The  first  term  represents  the  accumulation  in  volume  V  bounded  by  surface  S. 
The  second  term  is  the  net  mass  flux  through  the  boundary  surface.  Time 
averaging  the  above  equation  yields: 

/■f  [T//""]dv  +/[ff'fd,}"dSs0  (; 

V  T  ST 

Since  we  have  specified  a  flow  which  is  statistically  steady  in  time, 

l(fX  p*')’°  < 

I 

by  definition  and  this  result  can  be  generalized  to  any  function  <f>  as  well  as  p . 
The  mass  flux  term  can  be  used  to  define  a  new  integration  variable 

(i f7Pt i,)  t  ’  it  1 

The  sign  is  determined  by  the  product  U*  n  where "n  is  the  outer  normal  to  the 
boundary  surface  S.  The  continuity  equation  then  becomes 


J  dm  =  0 


For  a  mixing  problem  as  shown  in  Sketch  1 ,  the  equation  reduces  to 

•  »  ,  « 

=  m  (  + 

Note  that  continuity  has  not  defined  an  average  velocity. 
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The  momentum  integral  equation  is 

f  —  £UdV  +  J  (PUU-  r  -ndS  =0  (7) 

V  dt  S 

where  r  is  the  stress  tensor  defined  such  that  r  ."n  is  the  force  per  unit  area 
acting  externally  on  the  boundary  surface  S  and  the  sign  is  determined  by  the 
coordinate  system  positive  direction;  thus  the  normal  force  in  the  boundary  at 
station  1  is  positive  and  at  station  2  is  negative.  Normally  r  is  written  as  the 
sum  of  a  normal  stress  and  shear  stress 

*:-?!+  Ts  {8} 

Proceeding  to  time  average  the  momentum  equation,  the  first  term  of  Equation  7 
vanishes  if  we  define 

4>  *  IP^  )  (91 

fn  Equation  3,  From  this  point  we  will  consider  the  (streamwlse)  coordinate 
only,  uotiug  that  according  to  Sketch  1  the  outer  normal  n^  is  negative  and"n^  is 
positive.  Consider  the  momentum  flux. 

(~JVu  Udt*n)  *(  y  fu  (/5iU  £  )df  )  (10) 

T  x.  T  X, 


Equation  10  can  be  written  as  the  sum  of  two  terms,  the  mean  momentum 
flux  and  the  negative  of  the  normal  Reynolds  stress 

(fjfVuUn*.] i  *U^f  -<  V  "»  (l 


Since  by  definition 


>  |  /  -A  f  / 

ry  s  -  yj  <  >  u  d« 


The  time  average-  velocity  U  at  a  point  has  been  defined  by  the  momentum  flux. 
The  fluid  particle  has  only  three  degrees  of  freedom  at  a  given  point.  However, 
by  time  averaging  the  equation  of  motion,  there  are  also  effectively  3  internal 
degrees  of  freedom.  These  degrees  of  freedom,  if  excited,  contribute  to  the 
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mean  flow  momentum  and  energy.  The  mean  density  at  a  point  is  defined  as 
follows: 

—  Pu 
P  *  — • 

U 


thus  the  average  transverse  velocity  becomes 


V  = 


Vdt 


Averaging  Equation  7  yields 

"nl  dS  (13) 

X, 

Equation  13  can  now  be  integrated  using  the  boundary  condition  that  T>r  =  0  on  a 
solid  boundary.  There  may  also  be  area  changes  between  station  1  and  station  2 
on  which  the  normal  stresses  will  produce  either  a  thrust  or  drag  (T-D). 


f  f  4rfiP0U-f)n  dt  dS 

w$  T‘V 


]  .Ju  “»-/  [<V* 
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The  integrated  shear  stresses  will  be  expressed  by  a  friction  coefficient 
and  an  effective  L/  The  result  is 

m2  U2  +(P2  4-(?U)'u')2  A£=  m,  U,  +  (P  +(/>U)'u')1  A, 

+  mf  Uf  -f  (P  )f  Af 

+  (T-D)-4C(  ~  A,  (j  p  U2)  (14) 

The  velocity  in  this  equation  has  been  time  averaged  and  then  mass  flux 
averaged  and  thereby  the  mean  total  momentum  flux  has  been  defined.  The 
pressure  und  normal  Reynolds  stress  have  been  averaged  over  the  cross 
sectional  area  normal  to  the  Xj  coordinate. 

Note  that  turbulence  shows  up  in  the  momentum  equation  and  contributes 
to  the  stream  thrust.  Finally  the  energy  equation  in  integral  form  must  be 
time  averaged 
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fv  ~“<^«)<JV4-^/>UM-  ndS  =/  (U  - Ts)*n  <fS 

■*  s 

+jf  q  -7T  dS 

■»/  ~  dV  ( I5j 

\  df 

Now  by  time  averaging,  terms  1  and  5  of  Equation  15  will  vanish  and  the  viscous 
work  term  3  is  negligible.  The  fourth  term  is  the  net  heat  flux  which  is 
assumed  to  be  due  to  conduction  at  the  wall.  The  second  term  will  need 
expansion. 

By  definition  of  meau  and  fluctuating  quantities 


-  U  '  U  —  —  _  ll^t  |  *t  I  tli*t 

H=  h  +h  -P  4UU  fV  V'  4WW'4  (  ■  -T— _ ) 


then  the  integrand  of  term  2  becomes 


P  UH=  (^U)h  +(i°U)h'  +  (Aij  ■M/’UHlu' 

HP ff)vv’-M/>u  iww'  •» t pfi )( t5!r ) 

•ofc 

♦  +(PU)h'+(^u/ +{/=>U)'  u'U 

•HPU)  v'vfl/>U)  w'w  4</>U>'( — ) 

time  averaging  the  al>ove  expression  gives 
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/>UH 

=  PU  h  +■  PV 

I 

+  tJ  •  tt 

2T 

+  PS  {  ) 

( IS) 

ffl 

+  (/>U)'  h' 

3L 

+(/pu)*  (u*+v:w‘ ) 

The  turbulence  terms  in  the  equation  have  the  following  significance: 

I.  The  work  of  the  Reynolds  stresses  on  the  mean  flow 

II.  The  convection  erf  the  turbulence  kinetic  energy 
UI.  The  diffusion  of  specific  enthalpy 

IV.  Diffusion  of  the  turbulence  kinetic  energy 


The  equation  can  be  simplified  by  defining  the  turbulent  stress  tensor  to  be  the 
sura  of  a  hydrostatic  part  and  a  shear  stress.  The  turbulence  pressure  is 
defined  to  be  (for  incompressible  flow) 


PT  =  -j- 

Equation  19  can  be  simplified  to  read 


Pi)  H 


P  U  H 


where  H  and  q  are  defined  to  be 


and 


H 


h  + 


*»V 

u  •  u 
z 


+ 


PT  .  U,2+V>2  +W,Z 

P  2 


=  <  P  U  >v*  iP  U 


(20) 


(21) 


(22V 


(23) 
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Tbe  time  average  of  Equation  15  now  reads 


J  PVH  n  d5  =  j^U- l“T  *“s)n 


-  - - Js  ‘V  5 

/•  1%.  -te 

+  f  —  lq_  +q  )  •  n  dS 
•'S  T 

This  equation  is  now  integrated  over  tbe  control  volume  to  give  the  energy 
equation 


A  ,r  u^v»+w'f) 


2  2  2  P 


.  /  r  .U.-U,  pT,  U<2+v'2+W'2  x 

m,  V  h,  + - g  / 


+  rr<f  (hf  + 


,  %  .  U,2+V*+W*2 


2  r< 


+  Ch  pu,  (HW-HAW  )4-|-  a, 
From  Equation  22  a  turbulence  enthalpy  can  be  defined 

PT  .  U,2+V'2  +  W*2 

"T  =  “ .  >■ 


It  is  also  useful  at  this  point  to  introduce  the  concept  of  a  turbulence  temperature 
such  that  an  equation  of  state  exists  for  turbulence 


PT  =  p  rtt 

Thus  a  useful  expression  for  the  turbulence}  enthalpy  is 


{ 27) 


hT  *  -f-  «tt 


These  equations  reinforce  the  idea  that  turbulence  is  an  extra  internal 
degree  of  freedom  for  the  fluid.  For  thi*  degree  of  freedom  the  thermodynamic 
properties  correspond  to  those  of  a  monatomic  gas.  Turbulence  contributes  to 
the  thermodynamics  of  the  flow  but  with  a  separate  temperature  from  tbe  gas. 
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This  analysis  is  completely  general  at  this  point.  To  proceed  several 
assumptions  are  needed.  For  incompressible  flow  the  conservation  of  volume 
flow  rate  is  assumed.  First  let  us  calculate  the  turbulence  produced  by  a  baffle 
or  flame  holder  in  a  constant  area  pipe.  The  details  are  given  to  present  the 
method  of  solution  and  explore  this  basic  mechanism  of  producing  turbulence. 


jLL-L  tJJLL  Lf^lJJUJjLXMLJUJLLlXU 


— — > 


>  Turbulent 
State 


Sketch  2.  Control  Volume  for  Flow  Over  Flame  Holder 


The  conservation  equations  are: 
Incompressibility 


U*A2  =  U,  A, 


(29) 


Momentum 


m  U£*(P2  +PT  )At  *  m  U,  +  P,  A,  -D 


(30) 


Energy 


m  H-  -  rn  H, 


(31 ) 


where 


P_ 

P 


H  =  r  +  —  +  hT  4- 


u- 


(32) 


Assuming  that  no  dissipation  of  turbulent  kinetic  energy  has  occurred,  the 
internal  energy  is  constant.  Thus  with  Equation  29  the  energy  equation  yields 


hT  •  P,  -P, 


(33) 
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The  pressure  drop  can  be  found  from  the  momentum  equation.  Following 
Swithenbank  (Reference  2)  we  let 


®  "  2  P  ^n>a*  CQ  AB 


where 


u  .  -JL. 

umax  |  _  g 


Then  Equation  30  becomes  upon  substitution  of  Equation  29,  34,  35,  and  36 

P  :  "p  —  p"  — i_  piJ,2  — — - 

T2  **  Bg  I9"'  (|-Bg)« 

Thus  combining  Equation  37  and  33  v/e  find 


p  -  _L  pTr^fP-  ,f 

I  '  3  P  U>  8g  '  U,  *' 


or  we  can  write  the  turbulence  intensity  as 

-Ul  s  ( -L  -cP-bjl.\'4 

U  '  3  (l-Bg)*'  (39) 

We  can  also  determine  the  mean  static  pressure  urop 

3/1  — *  \  ®g 

p.  -p8  s  T  ^9pyi  )-■  --"  T  (40) 

832  (l-Bg)2 

These  formulas  Equations  39  and  40  are  not  the  same  as  Swithenbank  (Reference 
2)  presents  for  this  problem.  Comparing  specific  cases  we  find  for  B  =  0.25 
and  CD  =  0.4  that  Swithenbank  s  formulas  overestimate  the  turbulence  kinetic 
energy  by  nearly  a  factor  of  three.  This  analysis  shows  that  turbulence  is 
produced  by  a  velocity  difference  squared  while  it  might  be  thouijht  that 
turbulence  would  be  a  function  of  the  differences  of  kinetic  onerjjy. 


2.3  COMPARISON  OF  NORMAL  AND  COAXIAL  FUEL  INJECTION 


Next  we  will  look  at  two  mixing  problems  (Sketches  3  and  4j<. 


W UAtal m im dfUJuait  ii  hfAi^JUtall  IMkUXiJWM WU knh* 
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Sketch  3.  Control  Volume  for  Coaxial  Mixing 


u^I  ^  I u..Lj^ 


P  u  • 

ri  ui  , 


I  —  v 

t  \  ^ 


— .  I 

-_  I 

rrr  rr 
2 


A,  =A , 


Sketch  *1.  Control  Volume  for  Normal  Jet  Mixing 


These  are  two  practical  configurations  for  supersonic  mixing  as  well  as  the 
subsonic  incompressible  case.  It  is  of  interest  to  compare  these  two  to 
determine  the  effects  of  geometry  and  jet  velocity  on  the  turbulence  intensity 
and  total  pressure. 

2. 3, a.  Coaxial  Injection 

Assuming  the  fluids  are  incompressible  yields 

U2  A2  =  U,  A,  t-  TJf  Af  (41) 


Assuming  no  dissipation  of  turbulent  energy 


*2  =  Y,  f,  •+■  Yflf 


^  ,  -lr- 


(42) 


-Pf) 

i  r*  'rj  S7  L  "~r 

!'+-&)  {,  +  £} 


p,  -p.,  =p^+- 


-('+(■ 


I  Y  A< 
3  >  Af 


A 


)  j  (50) 
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2. 3.  b.  Normal  Injection 


For  the  problem  of  Sketch  4,  the  area  is  constant  and  there  is  no  fuel 
momentum  in  the  streamwise  direction.  The  solution  is 

_ £ 

U*  I  r  U,  2  U*  1/ 


where 


and 


U2=  U,  |+ 


f  A 


£  -  p 

2  rl 


'f 

(•+-£-) 


{(4- 


J_  A  ^ 


(51) 


(52) 


(53) 


Equation  51  can  also  be  written  in  terms  Of  the  known  inputs  and  solved  for 
turbulence  intensity 


For  the  normal  jet  there  is  always  a  pressure  drop. 

'P'|  -p»  *  \  *  f>,  [(  i+t)  ( l+  {-  -']  1  s5> 

Note  that  the  last  term  in  this  equation  is  constant  for  a  given  fuel  air  ratiOc 


2.3.c.  Total  Pressure  Losses 

In  the  high  speed  case  the  losses  or  gains  of  total  pressure  from  fuel 
injection  will  have  an  important  effect  on  the  total  performance  of  an  engine. 
It  is,  therefore,  useful  to  examine  the  total  pressure  losses  for  the  simpler 
incompressible  case. 


The  final  period  in  which  the  energy  decays  will  be  investigated  by  a 
similar  analysis  with  the  assumption  that  the  final  state  will  be  free  of 
turbulence. 
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4  >  :  ! 


Sketch  5.  Coatrd  Volume  for  Turbulence  Decay 


According  to  Sketch  5  which  represents  a  constant  area  section  in  which  the 
turbulence  decays  with  no  wall  friction  or  heat  transfer,  the  equations  are: 


Iucompressibility 


Momentum 


Uj  A  *  <^2  A 


n  U)  +  PjA  ~  m  U2  +  (P2+  P^) A  (57) 

Energy  (including  internal  energy  changes) 

P3  b,  „  P,  U*  ,  t 

«,  *-f  +-?  -  •«  +  -f  +  ~r  +  "t  <68> 

Solving  Equations  56  and  57  for  the  static  pressure,  we  find  a  pressure  recovery: 

P3  :  Tz  4-  PT  (59 ) 

From  Equation  58  the  dissipation  is  found  to  be 

«,  =e2  +  ~  RTt  (SO) 

Equation  59  shows  that  \*  is  the  sum  of  the  mean  static  pressure  and  the 
turbulence  pressure  that  should  be  considered  the  true  hydrostatic  normal  stress 
in  a  fluid. 

For  a  real  flow  the  processes  are  not  discrete  nor  one-dimensional,  and 
production,  convection,  and  dissipation  of.  turbulent  kinetic  energy  is  a  continuous 
process.  The  above  analysis  is  correct  for  state  5,  however,  and  ?s  useful  tor 
comparing  the  performance  of  rai  systems. 
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Let  us  now  compare  the  concurrent  and  normal  fuel  injectors  relative  to 
their  abilities  to  produce  turbulence  and  their  associated  pressure  losses  for  a 
given  iuel-to-air  ratio  as  the  jet  velocity  is  varied  while  holding  the  inlet  air 
conditions  fixed.  Pressure  losses  are  taken  to  be  the  difference  between  the 
inlet  total  pressure  and  the  outlet  static  and  dynamic  pressure  plus  turbulence 
pressure. 

While  these  results  are  essentially  only  applicable  to  incompressible  mixing, 
they  do  also  provide  some  insight  to  the  problem  of  supersonic  mixing  since  the 
velocity  ratio  will  be  of  the  order  of  unity  and  thus  density  changes  would  be 
small.  It  is  interesting  that  the  mixing  generates  turbulence  in  the  same  way  as 
a  baffle  by  velocity  difference  squared.  In  this  regard  the  normal  jet  is  slightly 
different  in  that  the  jet  kinetic  energy  is  not  controlled  this  way  and  it  is  eusier 
to  transform  jet  energy  into  turbulence. 

These  results  indicate  the  presence  of  influences  dus  to  geometry  and 
dynamic  conditions.  For  the  coaxial  jet  velocity,  ratios  between  0.5  and  2 
produce  relatively  little  turbulence.  The  normal  jet  produces  larger  turbulence 
intensities  with  smaller  variation  as  velocity  ratio  changes;  however,  no  thrust 
is  produced  by  normal  injection.  There  is  another  feature  of  the  normal  injector 
that  merits  note,  the  total  pressure  loss  is  constant  for  a  given  fuel-to-air  ratio; 
that  is,  it  is  independent  of  turbulence  intensity.  Thus  the  designer  cau  control 
mixing  rate  and  losses  for  this  type  of  injection.  It  further  points  out  that 
losses  are  determined  ultimately  by  the  boundary  conditions  and  not  by 
’’turbulence.*5  The  turbulence  is  the  intermediate  state  which  is  dependent  on 
the  total  dissipation  necessary  to  arrive  at  the  final  equilibrium  state  from  the 
initial  conditions.  Turbulence  should  be  considered  the  mechanism  of  tbG 
transition  processes  from  an  initial  unmixed  state  to  a  final  equilibrium  state. 
The  analysis  for  these  two  cases  will  be  extended  to  the  case  of  supersonic 
mixing  by  including  compressibility  effects.  This  extension  is  necessary  since 
’•flow  work?'  can  indeed  modify  the  results. 

2.4  COMPRESSIBLE  HOMOGENEOUS  TURBULENCE 

We  will  consider  the  ideal  situation  in  which  uo  transformation  of  turbulent 
kinetic  energy  into  thermal  energy  has  taken  place;  further  let  us  assume  an 
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absence  <>£  molecular  mixing  or  heat  conduction  which  is  consistent  with  the  first 
assumption.  These  conditions  represent  the  limiting  case  of  an  infinite  Reynolds 
number.  In  order  to  evaluate  the  integral  equations  of  motion,  we  will  have  to 
specify  the  spectrum  and  scale  of  turbulence.  However,  ii  appeal  a  that  this 
scale  5s  arbitrary  as  long  as  the  turbulence  is  homogeneous  and  no  gradients 
exist.  We  will  assume  that  a  definite  scale  for  the  "eddies”  exists  and  that  the 
turbulent  velocity  spectrum  of  these  is  single  valued  at  U'.  Consider  the 
following  sketch  of  a  set  of  cubical  cells  of  length  if  which  represent  the 
turbulent  eddy  flux. 


bt 


■£■  |f 

®  I 

I  i 

II  & 


f  f  I  f  o  I  f  j  f  a  !o  jf 


Sketch  6„  Compressible  Mixing  Model 


Some  of  these  eddies  are  fuel  eddies  and  some  are-  air  eddies.  For  the  i-th 


fuel  cell,  the  density  is 


V  p<  +  pr, 


The  velocity  can  be  written 


-U;  -  <U  +  U j )  T,  +  V',  tt  +  W'j  Ts 


where  by  assumption 


ff  #  If 

Uj  =  Vj  =  W}  :u 


Similarly  for  air  the  density  and  velocity  are 


^oj  s  Pa  +  Paj 


Uj  *  (U  +Uj  )«,  +•  Vj*  c“f  +  Wj  *, 
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Now  the  total  number  of  cells  which  cross  the  area  L  4  per  unit  time  are 

•  ■  -f 

where  total  number  is  the  sum  of  air  and  fuel  cells 


From  the  known  mixture  ratio 


which  can  be  reduced  to 


n  =  n0  -»-nf 


n? 

1  />fiL'S 

j=l  U  .  ± 

no  '  o 

2  P.  l'3 

Fi  0J 


J_  _L<l 

°  p< 


(67) 


The  mass  flux  per  unit  area  per  unit  time  is 

Pu  -  ~  J  PU4\  169) 

T 

or  considering  the  air  and  fuel  cells  separately  yields 

_  Of  _  no_ 

Pu  ■-  Y  £  /  (/>f +/>«  WU4-  Uj  >dt  +  --  2  J  ^0^^>U4-u'  )df  (70) 

*=«  *i  j=l  t, 

Now  we  assume  there  is  no  correlation  between  ,°'and  U*  then 

— rz-  _  _  nf  —  _  no 

P U  «  Pf  U  I  ».  +  />Q  U  I  tj  (71) 
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which  can  be  reduced  to 


and  similarly  for  the  air  eddies 


JL  J2l 

0  P, 


JL  Ji 

0  Pt 


i*4 

o  p 


Upon  substitution  into  Equation  72  the  mass  flux  becomes 

A.  -  P+f)Au 
A  P ♦*£) 


The  stream  thrust  can  be  treated  in  a  similar  manner.  The  expression  fi  r 
stream  thrust  is 


where 


Fr  -  (  />  U  U  )t  A  f  PA 


(  pu“u  )  =  —  jf  pu~u  df 


Considering  the  flux  of  eddies,  Equation  78  can  be  written 
-  r'f  _ 

(/>U  U  ~  t  J  +  +U,')(U+  u;  ) 

i=l 


!  na  -  _  _ 

♦f-  X  /  lPQ+PQ\  Hu+U-Xu+y!)  dt 

j=t  Tj  0  03  3  3 


which  becomes 


n»  n» 

r  .  r* 


(  />UU  ):  Pf  U*  I  t|  -*  />f  X  u;  ». 

i=i  i=l 


_ ,  "o  _  "n  . 

+  po  u  I*,  fi  ¥  V*  u'ffi 


19 


AFAPL-TR-71-18 


Now  since  U^  =  U  ^  this  equation  reduces  to 


+  PQu  —  +  p0  u  -f 

Now  Equation  81  becomes  by  substitution  of  Equations  74  and  75 

.2 


(  PU U  )  = 


^Ud  +  4-)  _  A.U*  (l  +  — ) 


-=?—  u  + 


(,+  X.  JjL\ 
\  a  pf  i 


(,-f-f  A) 
a  Pf  1 


The  stream  thrust  Equation  (77)  becomes 


(8! ) 


(82) 


<V  (  P  +PT)  A  +  m  U  (83) 

Again  the  turbulence  contributes  to  the  stream  thrust.  Finally  the  energy  flux 
will  be  obtained  by  the  same  method. 


The  energy’  equation  requires  au  average  enthalpy  flux 


PU  H  z  yIt  pV  H 


dt 


(84) 


where 


_  ,2  ,2  ,2 

,  u2  —  .  0  +v  +w 

H  =  h  -i  n  -b  +  U  U*  +•  - g - 


(85) 


Substituting  Equation  85  into  Equation  84  and  separating  the  air  eddies  from  the 
fuel  eddies  yields 


(PUH)X=  Y  f  /( pf+/>;.  hufu!  )(hf+  hu  f  ^+uu.'+ f-  u!* )  at 

i=i  t i 

m  9  —~2 

+  4'  £  Ij  VVtIT+ui  *  ~+ua+|-uf)dt  (86) 

J1!  ‘  C  \  C.  \ 
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2.5  SOLUTION  OF  THE  COMPRESSIBLE  EQUATIONS  OF  MOTION  FOR 
COAXIAL  AND  NORMAL  INJECTION 

We  will  consider  gaseous  injection  for  the  coaxial  and  normal  fuel  injector. 
Then  let  us  further  assume  the  fuel  is  a  diatonic  gas.  Thus 

ra  z  rt  =l-4  06) 

We  further  assume  that  the  fuel  is  injected  at  the  same  static  pressure  as  the 
air  without  creating  shock  waves 

P,f  =  (97) 

Let  us  introduce  the  following  nondimensional  variables 


Then  the  equations  become  for  a  unit  mass  flow  of  air 

(I+17)  u‘‘  P  ~r  =A2  {105) 

PAf(l  +  |!{^  -l)U  (HU'*)  =FZ  (i06) 

i\+V)  jrVy 

- : —  P  +(0.-t)U  (1  +  5U/  )-  9  -0  007) 

(H-fl 

If  we  introduce  a  new  variable, 

Xr  =  P  i  108) 

We  can  reduce  Equations  102,  103,  and  (104)  to  a  single  equation  in  X 

A12  X-  t-AB  X10**,  X5  *A0  =0  (109) 

where 


ti+77 )  ,  9.  —  \ 

<!+-*•  i  ^  F,->  ' 

(  110) 

Aio  ■  -  Bt 

(til  ) 

(6,-n  n-i-77)  f2 

A5  -  5  * . 

(  F,  -  1)  (i  +-— )  A  2 

(112  ) 

A0  =  -4(0,  -1)  (I+tj)2  a| 

(113! 

f,  =  1+yMf 

(114) 

e,  *  h-  «f 

(115) 

A  =  |+«  *  C7?*  T°f)  (-^-)  °' 

A*  1  *  «  T0s  '  'P0f  '  Of 

(116)' 

*See  list  of  symbols  for  definition  of  D,,  D,, 

X  * 
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These  conditions  result  in  a  velocity  ratio 


2.24 


Comparing  this  result  \%  5th  that  obtained  for  incompressible  flow  (see  Figure  2), 
we  find  a  significant  effect  of  compressibility  which  reduces  the  amoum  of  energy 
transformed  into  turbulence.  However,  it  appears  that  there  is  still  a  sufficient 
supply  of  energy  for  turbulent  mixing  to  predominate  over  molecular  diffusion. 
This  reduction  in  the  maximum  energy  of  turbulence  can  be  one  reason  to  explain 
the  observed  slower  mixing  rates  in  compressible  flows. 
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The  effect  of  the  parameters  M.,  -5^,  and  —  are  shown  in  the  Figure  3. 

Pof  1  P0l  a  XT. 

The  effects  of  Mj  and  -jj—  are  minor.  Thus  we  find  thatW-  is  the  most 

01  ua 

important  dynamic  parameter  for  the  compressible  case;  however,  it  is 

controlled  by  the  temperature  and  molecular  weight  of  the  fuel  entering  the 

combustor.  Thus  we  find  it  is  important  to  simulate  correctly  the  fuel  total 

temperatures  in  combustor  development  work  with  boilerplate  engine  hardware 

in  order  to  obtain  the  correct  mixing  rates  and  heat  release  rates  which  may 

interact  with  each  other. 


2.6  FROZEN  FLOW 


It  is  of  fundamental  interest  to  calculate  the  thrust  loss  in  a  supersonic  nozzle 
due  to  a  frozen  pattern  of  turbulence.  We  will  assume  that  the  gases  mix  on  a 
molecular  scale  so  that  thermal  equilibrium  prevails  Hat  let  the  turbulence 
kinetic  energy  be  frozen.  Further  we  will  assume  that  this  molecular  mixing 
takes  place  in  a  constant  area  section;  heat  is  added  downstream  and  then 
expanded  in  the  nozzle.  The  process  is  depicted  in  the  following  sketch. 


Sketch  7.  Control  Volume  for  Turbulent  Combistion  and  Expansion 


The  energy  equation  can  be  written  as 


—  —  y  2 

H  =  h  +  hT+  -y 


(120) 
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Since  for  a  perfect  gas  y  -  1.4,  this  equation  can  be  reduced  to 

T0  =  T  (ltyW*H-^r  Tt 


Let  us  define  an  effective  total  temperature 


tot  =  to  — y  tT 


>t  *  ?  i,+2T-uZ) 


(123) 


If  the  gas  is  brought  to  rest  isentropically  while  the  turbulent  motion  remains 
the  same,  we  have 

y 

*0T  =P  («+  (124) 


To  simplify  the  formulation  of  the  problem,  let  us  assume  that  the  properties 
of  the  comtcf.tion  products  at  t  tation  3  do  not  differ  from  those  of  air.  Then 
conservation  of  mass  and  momentum  yield  the  following  equations  for  Mach 
number  at  station  3 


M«  (■+  f  ('+4xt°tJ4 


where 


3  '  A. 


(IH.  -IL  +  4-  y  Mj  )”* 


r°3  s  (T°,  ~cT)  /(,  +  f  ) 


(125) 


(126) 


(127) 


Now  there  is  also  no  loss  of  generality  In  assuming  that  the  expansion  is 
iseotropic  when  the  turbulence  is  frozen.  The  equation  for  Mg  is 
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(74-1)  , 

2(y-D  (i^j-  )(T°Ti) 

'V  {xp 


Now  since  the  stream  thrust  is 


-  /  TT 

F«  *  P.  (  I  +  -~- 

we  can  write  for  specific  impulse 


( 


(  128) 


(129) 


(130) 


This  analysis  allows  us  to  compare  "equilibrium”  performance  and  "frozen" 
performance.  Frozen  expansion  losses  are  depicted  in  Figure  4.  When  the 
results  are  plotted  in  this  manner,  there  is  no  difference  between  coaxial  and 
normal  injection.  A  small  effect  was  observed  when  heat  was  added;  the 
results  presented  are  for  the  limiting  case  of  constant  total  temperature. 


These  losses  seam  small  enough  for  the  large  turbulence  temperatures 
involved.  First,  the  turbulence  pressure  contributes  to  the  stream  thrust  and 
contributes  to  the  flow  work  in  the  nozzle;  second,  a  large  part  of  this  energy 
is  unavailable  in  any  case  since  it  will  generate  entropy  when  dissipated.  In 
fact,  the  total  pressure  decreases  when  the  turbulence  is  relaxed  in  the  constant 
area  process.  These  performance  losses  are  nearly  one  order  of  magnitude 
smaller  than  presented  in  Reference  1.  Thus  this  fundamental  .limitation  to 
performance  of  supersonic  combustion  ramjets  is  not  as  severely  restrictive 
as  previously  thought. 


2.7  RELATIVE  PERFORMANCE 

Finally,  let  us  compare  the  effects  of  'uel  injection  mode  on  ">erfcrmance 
of  the  engine  when  the  air  and  fuel  initial  conditions  are  identical  sr.cept  ior 
geometry.  For  the  coaxial  mode  the  fuel  momentum  contributes  directly  to  the 
engine  thrust.  The  normal  jet  contribute?  no  momentum  and  the  air  must 
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trausfer  momentum  to  the  fuel  stream;  this  behavior  results  in  an  overall  Mach 
number  reduction  and  total  pressure  loss.  However,  the  heat  audition  losses 
will  be  reduced  at  this  higher  pressure-  but  lower  Mach  number  condition.  A 
larger  expansion  ratio  is  also  available  from  the  combustor  to  nozzle  exit  since 
there  is  no  area  change  in  the  mixing  region  as  in  the  coaxial  mode.  The 
relative  effects  depend  on  the  fuel  total  temperature,  total  pressure,  and 
molecular  weight. 

Results  presented  in  Figures  5  and  G  show  that  very  large  mixing  losses 
occur  for  the  normal  jet  as  compared  to  the  coaxial  jet.  However,  total 
pressure  losses  with  heat  addition  are  comparable  and  the  loss  of  engine 
performance  is  small, 


2.8  SUMMARY  OF  DESIGN  CONSIDERATIONS 

This  analysis  comparing  the  characteristics  of  supersonic  mixing  has 
yielded  the  following  results  which  are  important  design  considerations. 

1.  Fuel  injection  mode  is  the  biggest  single  factor  affecting  turbulence 
intensity,  normal  jets  produce  much  more  turbulent  energy  than  coaxial  jets. 
Thus  normal  jets  should  have  faster  mixing  rates. 

2.  Normal  jets  produce  nearly  constant  turbulence  intensity  for  a  giveu 
fuel-air  ratio  independent  of  other  initial  conditions  while  the-  energy  produced 
by  coaxial  jets  varies  widely  with  jet  velocity  ratio.  Thus  normal  jets  should 
have  more  constant  mixing  lengths  as  flight  speed  varies. 

3.  Frozen  turbulence  performance  losses  are  not  large  but  do  increase 
with  flight  speed.  Thus  large  turbulence  intensities  can  be  generated  to 
achieve  rapid  mixing  at  all  but  extreme  hypersonic  speed. 

4.  Normal  jets  do  not  produce  large  performance  losses  compared  with 
coaxial  mixing;  however,  a  constant  area  burner  will  choke  at  a  lower  flight 
speed  when  normal  injection  is  used.  Thus  coaxial  jets  may  be  very 
advantageous  for  low  flight  velocity. 


It  is  concluded  that  normal  jet  fuel  injection  that  maximizes  turbulence 
intensity  should  be  developed  in  order  to  minimize  and  stabilize  mixing  length. 
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SECTION  HI 

TURBULENT  FIELD  EQUATIONS  FOR  COAXIAL  MIXING 


3. 1  COMPRESSIBLE  MIXING 


In  the  mixing  field,  turbulence  is  produced  and  dissipated  continuously. 

There  is  a  strong  coupling  between  the  mean  flow  and  turbulence.  The  mean 
flow  determines  the  scale  of  turbulence.  The  scale  and  the  mean  profile 
determine  both  the  rate  of  production  and  the  dissipation  rate  of  the  turbulent 
kinetic  energy.  The  coupling  is  complete  since  this  turbulence  intensity 
directly  determines  the  mean  velocity  field  through  the  Reynolds  stresses.,  in 
the  approach  that  follows  a  unique  relation  between  turbulence  intensity 
(turbulent  kinetic  energy)  and  "eddy*1  viscosity  is  assumed  to  exist.  Further  a 
unique  relation  between  the  intensity  and  scale  of  turbulence  and  the  dissipation 
rate  is  also  assumed  to  exist. 

In  this  regard  the  turbulence  must  be  fully  developed  such  that  an  equilibrium 
spectrum  function  exists  or  nearly  exists.  Thus  these  assumptions  and 
formulation  which  follow  are  valid  for  this  "equilibrium’’  turbulent  flow  field. 
Thus  the  question  of  intermittency  is  raised.  No  direct  account  of  intermittency 
is  taken  in  the  present  formulation  since*  where  intermittency  is  large,  the 
turbulent  transport  is  small  and  the  overall  effect  on  the  flow  field  is  small. 

Hear,  transfer,  diffusion,  and  momentum  transport  are  all  assumed  to  be 
similar.  It  is  further  assumed  that  the  eddy  viscosity,  Prandtl  number 
relation,  etc. ,  iound  in  any  turbulent  flow  will  be  valid  in  all  turbulent  flows  as 
long  as  the  above  restrictions  are  true. 

3.2  TURBULENT  STRESS  TENSOR 

The  Reynolds  stresses  are  defined  by  the  following  expression 


ty  *  -  (PU)‘  u' 


AFAPL-TR-7 1-18 


which  for  cylindrical  coordinates  has  the  following  components; 

rrr  s  -  1/>U)'UT  (  132) 


t,,  - 

'from  dynamic" 

(133) 

T'Q  --W'O'  *  TQf 

"equilibrium" 

034) 

TBi  ~W,V  "T,S 

(135) 

rQ9  -</>V)V 

(136) 

Ttz  w< 

(137) 

Let  u&  introduce  the  eddy  viscosity  so  that  we  can  write  the  turbulent 

2 

stresses  as  the  surn  of  a  normal  "hydrostatic”  part  and  a  shear  stress 

rT  =  -  PT  I  +  fij  def  U  (138) 

where 

*  P*  139) 

such  that  c  is  a  kinematic  eddy  viscosity. 

The  expression  for  turbulence  pressure  is 

PT  *  PyK  (140) 

where 

K  =  Y*  (  Trr  +  XQQ  +  TJJt  )  /  P  (!4I) 

is  the  turbulent  kinetic  energy. 

■■  .—>■■■  ■■),■■■■■  ■»■■■■■ 

2  The  deformation  tensor  defU  isfjjrmed  from  the  V  operator;  defU  :fo  +  { VU)* 
wbera  {$$)*  is  transpose  of  (VU). 
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Then  the  components  can  be  written  as: 


r„  -  2“[ 

x  .  au  ] 

"  3  ^  *  dr  J 

(142) 

roe  ‘  2p{ 

3  *  €  '  r  dQ  +  r  )  j 

(143) 

T:i  ‘  Zp\ 

k  dH  i 

■t  4  •  Tr  J 

(144) 

Trz  z  ~P*  ( 

dV  ,  d  U  , 

‘  df  dz  1 

(145) 

Tr6  s  P*  ( 

r  d  (  V>  i  au  , 

dr  [  '  >  '  66  1 

(146) 

vdz  ,  P*  1 

>  Jly.  +  jl  ii\ 

1  dz  r  69  1 

(147) 

It  can  be  argued  that  the  eddy  viscosity  should  be  a  function  of  both  the  mean 
flow  scale  and  the  turbulence  intensity.  Obviously  the  eddy  sizes  are  much 
larger  in  ocean  currents  than  in  laboratory  flows  and  thus  the  cross  velocity 
products  should  be  correlated  for  much  longer  periods  of  time,  while,  if  the 
turbulence  intensity  vanishes,  the  cross  correlations  also  vanish.  It  will  be 
necessary  to  determine  this  relationship  experimentally. 


3.3  DYNAMIC  EQUATIONS  FOR  MEAN  MOTION 


The  differential  equations  for  the  time  mean  flow  are  well  known.  We 
present  them  here  without  derivation  in  cylindrical  coordinates  for  flow  with 
axial  symmetry  in  the  usual  boundary’  layer  approximation. 


Continuity 


T~ilrlr  ^  1  4  7T  * =0 


(148) 


Diffusion 


J-j-lrpUY;  1+^ipwY;  +<D  j  (149) 
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where 


iP  U)'  Y,' 

I-* 


dr 


(150) 


Radial  Momentum 


T  4-  ( U)  U  =  P_ 


(151) 


Streamwise  Momentum 


7- (r  PU  W)  +  3^- (/>W  W  )  =  -  +  7- )  I7-) 


(152) 


where  we  have  assumed 


PT  =  (  />U)  u'  =  (  P  W)'  W 


153) 


Total  Enthalpy 


7-§7  <" pu  H,+  * 1  *  T-  ■£ 


dr 


dY; 

\ 

dr 

1 

.  1 

d( 

+  r 

dr  ' 

-  1 

dr 


(154) 


where 


(^U)V  =  Cp  (PU)'t'  +  I  hj  (/CU)'  YjJ 


^TTT 


(155) 


With  the  definitions  of  the  turbulent  diffusion  coefficients  for  heat  and  mass 
transfer,  v/e  obtain 


(PU)'  h'  =  -P CP  €T~  -  I  hj 
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Then  the  energy  equation  can  be  written  as: 


i  d 
'  dr 

( r  PU  H  )  +  (  pW  h  1  =  -7* 

dz  r 

it- 

/^«\  dH  1 
'Pf‘  dr  i 

4-j-f  r 

r  dr  L 

-($•1 

+4  —  f 

r  (-£?-) 

.  »  WZ/2 
KP  -1 )  — - 

r  dr  L 

\  Pr  t 

^  r  dr 

+  ±  ±-  rr  /JjS.v3f!L -A 

r  dr  L  '  pr  ^5  «T  )  dr  J 


dhi 


Turbulent  Energy  Field 


The  turbulent  kinetic  energy  equation  becomes 


UK  >  + -Jj  CW  K  1  .  M  (-^L)  -  PT 


dk 


where 


dk  =  d  +•  U'-V  P# 


(167) 


The  terms  are  explained  as  follows: 

d  dissipation  of  turbulent  kinetic  energy 

U.#V  P*  acoustic  coupling;  this  term  is  assumed  to  be  a  small  loss  unless 
an  interaction  with  a  sound  source  is  introduced.  Such  effects 
have  been  observed  and  can  be  large. 


■J 

i 

i 


An  assumption  has  been  made  regarding  the  role  of  viscosity  on  the 
generation  of  turbulence.  We  have  assumed  that  this  action  on  turbulence  is 
purely  dissipative.  Thus  the  viscous  production  term  in  the  equation  is  replaced 
by  the  turbulence  dissipation  rale 

=  -d  { 1681 

and  in  the  total  enthalpy  equation  the  turbulent  viscous  woik  becomes 

7  *  (  U'-  r's  1:0 

34 
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The  arguments  for  this  assumption  are  discussed  below. 


Dissipation 


The  dissipation  of  turbulent  kinetic  energy  *n  isotropic  incompressible 
turbulence  has  been  shown  to  obey  the  law  {Reference  8), 


U 


/» 


2  ^  Ld 


( iroj 


where 


L~  =  the  dissipation  scale  which  is  approximately  equal  to  the  integral  scale 
(or  scale  of  large  "eddies'*} 


This  law  determined  by  dimensional  reasoning  and  verified  by  experiment  can  be 
explained  by  existence  of  an  equilibrium  turbulence  vekcity  distribution  function. 
The  dissipation  occurs  at  the  smallest  scale  while  the  energy  of  the  smallest 
scale  is  directly  related  through  the  distribution  function  to  the  large  scale 
eddies  which  contain  most  of  the  energy.  This  dynamic  equilibrium  determined 
by  the  energy  transfer  process  from  one  scale  to  another  is  essentially  viscosity 
independent  since  the  energy  flows  from  large  scale  to  small  scale  and  is  finally 
terminated  by  viscous  dissipation  into  heat  at  the  end  of  this  process.  Since 
data  are  available  for  incompressible  flow  only,  the  dependence  of  dissipation 
on  turbulence  intensity  will  be  determined  experimentally. 


Static  Enthalpy 


The  equation  for  static  enthalpy  is  somewhat  simpler  to  use  for  solutions  by 
numerical  integration  and  can  be  written 


— —  d  h 
+■  i&W  ~ - -  W 


dz 


dz  * 


1  Hfftc 


i  o  r 
r  dr  V 


Pr 


dJ 
P  dr 


) 


,  1  d  t 
+  r  dr  'r  Sc 


rt  iXL\ 

i  IP 


(171) 
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where 


.  ...  -/aw  \* 

*•  s  dk  +  M  l—) 


(172) 


The  form  of  these  equations  makes  them  applicable  to  laminar  as  well  as 
turbulent  flows  when  K-M).  Since  the  boundary  layer  equations  are  parabolic, 
we  must  specify  the  pressure  gradient  in  order  to  obtain  a  solution. 


3.4  TRANSFORMATION  OF  EQUATIONS 


The  equations  will  be  solved  in  the  Von  Misas  plane  by  using  the 
transformation 


Z  =  2 


Mr  s  />Wr 


ffz  =  -  PU  r 


The  equations  of  conservation  then  become 


Momentum 


=  W  dZ  dr  dti/V 


dZ  PW  dZ  \Jr  df'  df 


aw 

°  4  dr  ) 


Mean  flow  energy 


1  a  P  ,  i  a  t cp 0  dT 


-  1  V  ^  ^  /  ~  K  V  »  \ 

az  p  az  \y  dy '  Pr  adr  * 


+J fL(-2_  y  ^  -^Yi  >  ,  la, 

ifr  dty  '  Sc  }  •  a^  '  p  W 


Turbulence  energy 


+  j.  a  t  o  ^<k  ax . 
^  a^ '  a^  * 


c  ,  aw 


pt  aw 


+  f  iTf)  ~  Ji 


(175) 
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Diffusion 


d  Yj  _  d  /  o  d  Yj  \ 

dZ  dt  '  sc  dt  ' 


where 


(178) 


(179) 


3.5  NUMERICAL  TREATMENT  OF  DIFFERENTIAL  EQUATIONS 


f’Ketch  8.  Definition  of  Differences 


This  finite  difference  formulation  for  the  solution  of  the  turbulent  equations  of 
motion  is  a  minor  modification  of  the  work  of  Edelman  and  Fortune  (Reference  9), 
The  definition  of  the  f  .nite  difference  analogues,  for  any  variable  <f>  of  the 
differential  operators  in  the  equations  of  motion  are  defined  according  to  Sketch  S, 
as  follows,  ’’’he  difference  between  this  formulation  and  that  of  Edelman  is  the 
inclusion  of  L«c  v.  ria  ion  of  the  eddy  viscosity  across  the  mixing  layer  in 
Equation  182.  The  first  derivative  in  the  Z  direction  is 


.ii.) 

dl  V.*l,nr» 


S7 
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Diffusion 


Y|n+I,m  'Y'n,m  +  ^  Sc -2Y'n,m  +  Y'n,m-{) 

The  static  enthalpy  equation  written  in  terms  of  static  temperature  becomes 


/  (  dP  \ 

AZ 

/  dk  \ 

\  vaz  'nfl 

P„,  m  + 

'  />W'n 

AZ  ,  £po  x  /_  ^  ^ 

“2Tn,m  +Tn,m-lj 

+  f&<F‘  “  Y'n,<"-l)  (Tn.'f.+|  "Vir-l)]}/ 


(CD  ’ 


P  P 


•  n,m 


(186) 


where 


*  /  dP  \  .  /dPo>  (  /  . „  \  2  .  Py  . ,  dfc  \ 

/J,  m  '  dZ  'n+|  '^2  A»+l  flips'  p«  'n,m  3  v1  4  Klz'n  +  I, 


3,6  EQUATIONS  OF  MOTION  ON  THE  AXIS 

These  equations  have  a  singular  point  on  the  axis  as  r  -~*Q  or  \fr—*0  and 
their  limiting  form  must  be  investigated.  This  limiting  form  is  also  useful  for 
evaluating  experimental  data  to  deduce  or  determine  inversely  the  turbulent 
transport  coefficients. 

Streamwise  momentum 

^s2,e4.i!i.  (.88) 

r  di  ^*dr  dz 
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Static  enthalpy 


P * 


dh 
d  z 


Turbulence 


Diffusion 
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The  equations  in  the  Von  Mises  Plane  are: 
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Finally  the  difference  equations  ar .. 
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Energy 


w  =V'  +  x.(^^"AZ 


Turbulence  kinetic  energy' 


<n  +  ».'  *  Kn,r  AZ  (  r>w) 


pW/ntJ 


Diffusion 
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3.7  THE  INVERSE  TRANSFORMATION 


The  physical  coordinates  (r,  z)  can  be  found  by  the  equations: 

z  =  2 

~o 

Equation  200  expressed  in  difference  form  can  be  written 


rnf  i,m  =[fn+i,fn-i+A^  ("TI 


m 


Tn-t 


<^W,n4-t,m  lPwin+,t  m-‘/ 


3.8  SUMMARY  OF  THEORETICAL  CONSIDERATIONS 


In  Section  n  it  was  shown  that  turbulence  has  an  equation  of  state 


Py  =  P  RTy 


( 196) 


(197) 


(198) 


(199) 

(200) 


(201) 


(202) 


and  that  associated  with  this  turbulent  state  were  the  state  functions  of  turbulent 
energy  and  turbulent  enthalpy.  In  the  present  section,  it  was  assumed  that  the 
transport  properties  of  the  turbulent  flow  are  also  functions  of  the  turbulent  state. 
Thus,  it  is  assumed  that  an  eddy  viscosity  exists;  that  io,  that  eddy  viscosity  is 
functionally  related  to  the  equation  of  state.  This  assumption  implies  (hat  an 
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equilibrium  turbulent  velocity  spectrum  exists  such  that  shear  produces  only  a 
minor  departure  from  this  equilibrium  structure. 

The  dissipation  into  thermal  energy  f^r  turbulent  shear  flow  must  also  be 
related  Co  the  turbulent  state  under  these  assumptions.  Thus,  the  experimental 
determination  of  the  relations  between  eddy  viscosity  and  turbulent  energy  and 
between  energy  dissipation  and  turbulent  energy  are  required  to  verify  the 
present  theoretical  model. 


it 


^'gcZSa&g 


e 


AFAPL-TR-71-18 


SECTION  IV 

DESCRIPTION  OF  EXPERIMENTS 


4.2  SUPERSONIC  FREE  JET 

The  objective  of  these  experiments  was  to  obtain  basic  data  on  mixing  rates, 
dissipation,  and  turbulence  intensity  in  a  compressible  flow.  The  jet  Mach 
number  and  total  temperature  were  varied  independently  by  use  of  interchangeable 
nozzles  fitted  to  a  hydrogen  fueled  burner.  The  total  temperature  was  restricted 
to  1500°Il  in  order  that  the  combustion  products  would  not  deviate  appreciably 
from  calorically  perfect  air. 

4.1.  a.  Nozzles 


The  nozzle  internal  contours  were  computed  by  the  method  of  characteristics 
for  a  y  ~  1.4.  The  nominal  design  Mach  numbers  of  1.0,  1.5,  2.0,  and  2.5 
were  chosen  based  on  limitations  of  the  available  air  supply.  The  nozzle  axit 
diameters  were  all  equal  to  0. 1  foot .  Exit  Mach  number  profiles  are  shown  in 
Figure  7a  through  7d.  Typical  total  temperature  profiles  at  the  nozzle  exit  are 
shown  in  Figure  8a  through  8d. 

4, 1.  b.  Instrumentation 

The  free  jet  combustor  was  instrumented  for  pressure  and  temperature. 

The  jet  was  probed  for  impact  pressure,  static  pressure,  and  total  temperature. 
Consolidated  Electrodynamics  Corporation  pressure  transducers  with  bridge 
balancing  networks  were  used  to  measure  pressure  for  temperature  measure¬ 
ments.  Iron  constantan  thermocouples  were  used  with  Brown  Instrument 
vertical  scale  potentiometers  and  a  room-temperature  reference  junction. 

Probe  position  was  determined  by  calibrated  voltage  drops  from  battery-powered 
variable  resistors  mounted  on  the  probe  positioning  cart.  All  data  were 
recorded  on  magnetic  tape  with  a  digital  data  acquisition  system  SRL-200  scanning 
at  a  rate  of  2G0  channels/  second.  Schlieren  photographs  were  taken  with  a 
spark  light  source  of  0. 2-microsecond  time  constant.  Total  pressure  and 
temperature  surveys  of  the  jot  were  obtained  with  a  combination  probe  shown 
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in  Figure  9a.  This  probe  was  not  explicitly  compensated  for  radiation  or 
conduction  losses.  Various  size  thermocouple  wires  were  tested  until 
temperature  readings  at  identical  conditions  did  not  change.  A  5-mil  thermo¬ 
couple  wire  was  found  to  give  good  results.  Radiation  compensation  was 
provided  for  by  surrounding  the  thermocouple  with  a  low  conductivity  material 
(LAVA)  which  minimized  the  radiation  losses  from  the  thermocouple.  Probe 
bleed  was  maximized  to  provide  as  rapid  a  response  to  changing  conditions  while 
not  degrading  the  total  pressure  measurement.  A  simple  impact  probe  was  used 
to  verify  the  pressure  measurement.  Temperature  measurements  were  also 
made  with  a  conical  recovery  temperature  probe,  developed  at  the  Naval  Orduance 
Laboratory,  White  Oak,  Maryland  (Figure  9b).  These  measurements  were  in 
agreement  with  those  made  with  the  combination  probe  (Figure  10). 

Static  pressure  surveys  were  made  with  the  probe  shown  in  Figure  9c. 

Three  static  taps  were  positioned  120°  apart  at  ihe  same  axial  location.  It  was 
hoped  this  spacing  would  allow  the  true  static  pressure  to  be  measured 
uninfluenced  by  the  cross  flows  induced  by  the  turbulent  velocity  fluctuations. 


4.  l.c.  Procedure 


Jet  operating  conditions  were  established  and  allowed  tc  stabilize  as  near 
the  design  operating  conditions  for  the  nozzle  as  possible.  Then  radial  impact 
pressure  and  total  temperature  surveys  were  made  at  three  axial  positions 
followed  by  a  centerline  survey  of  impact  pressure  and  total  pressure.  Next, 
the  static  pressure  probe  was  positioned  on  the  jet  centerline  and  a  centerline 
survey  of  static  pressure  was  recorded.  Since  the  data  system  was  of  the 
analog  to  digital  converter  design  which  takes  essentially  instantaneous  readings, 
10  scans  wore  taken  at  each  position  so  that  most  of  the  noise  could  be  averaged 
out  of  the  pressure  measurements. 

4.  l.d.  Accuracy  of  the  Measurements 


Pressure  transducers  used  to  measure  total  pressure  were  rated  at  1/  2% 
of  full  scale  accuracy  and  were  changed,  to  give  the  largest  possible  voltage 
reading,  for  each  Mach  number  nozzle.  Thus  good  accuracy  was  obtained  near 
the  early  part  of  the  jet  decay  but  it  decreased  as  the  decay  increased.  The 
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static  pressure  transducer  had  a  0-25  PSIA  range  and  the  accuracy  should  hove 
been  constant  in  the  whole  flow  field.  No  attempt  was  made  to  determine  any 
probe  effects  on  the  accuracy  of  the  static  pressure  probe  other  than  tc  introduce 
small  amounts  of  deflection  in  the  proble  and  compare  results  from  different 
experiments.  No  trends  were  observed  that  could  not  be  attributed  to  run 
variation  in  jet  conditions.  Another  source  of  error  waa  that  of  bridge 
unbalance  which  had  various  causes.  This  error  had  tho  effect  of  shifting  the 
calibration  a  constant  amount  over  the  entire  range.  This  effect  was  corrected 
by  taking  a  reading  with  no  flow  and  comparing  it  with  the  1  atmosphere 
calibration  values  for  each  transducer. 

As  an  overall  system  check,  mass  flows  were  computed  from  measured 
chambers  conditions  and  were  compared  with  iniegraded  radial  profile  surveys 
taken  at  the  nozzle  exit.  The  results  are  shown  in  Figure  11.  The  agreement 
is  good;  the  maximum  deviation  being  about  15%.  Shown  in  Figure  12  is  the 
comparison  of  integrated  with  calculated  exit  stream  momentum;  the  agreement 
is  again  quite  good  with  the  maximum  deviation  being  about  6%.  Integrated 
excess  enthalpy  and  stream  thrust  from  various  experiments  at  downstream 
stations  are  shown  in  Figures  13  and  14,  These  integrals  are  theoretically 
constants  of  motion;  however,  some  noticeable  deviations  from  constancy  are 
seen.  The  large  deviations  are  believed  to  be  caused  by  the  reduced  accuracy 
of  the  probe  position  measurement  at  large  distances  from  the  centerline  where 
the  signals  are  low  but  make  significant  contributions  to  the  integrated  values. 

4.2  COAXIAL  SUPERSONIC  MIXING 

The  objective  of  these  experiments  was  to  determine  the  structure  of  the 
mixing  field.  A  hot  film  anemometer  was  used  to  determine  the  extent  of  the 
turbulence  field.  Helium  was  injected  coaxially  from  a  Mach  4  conical  nozzle 
into  a  Mach  3  air  stream.  The  gas  injection  nozzle  was  attached  to  a  half-inch 
diameter  steel  tubing  which  extended  through  the  facility  nozzle  into  tfce  test 
section.  The  long  length  of  the  steal  tubing  gave  rise  to  a  thick  turbulent 
boundary  layer  on  the  outside  of  the  tubing.  Since  it  was  impossible  to  remove 
the  boundary  layer,  its  effect  on  mixing  was  not  determined  experimentally. 

The  wind  tunnel  had  a  free  jet  test  section  and  a  Mach  3  nozzle  of  rectangular 
cross  section,  4  in.  x  7  in.  The  usable  test  section  length  was  about  12  in. 
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Mixing  of  the  gas  of  this  jet  with  the  surrounding  air  was  found  to  be  extremely 
slow.  The  jet  was  operated  so  that  the  jet  exit  pressure  was  equal  to  the  air 
static  pressure. 

4. 2.  a.  Instrumentation 

Instrumentation  for  this  experiment  consisted  of  pressure  transducers  and 
thermocouples  to  measure  upstream  reservoir  pressure  and  temperature  of  both 
air  atxl  helium.  Impact  pressure,  static  pressure,  and  output  of  a  conical  hot 
film  anemometer  probe  were  also  measured.  The  hot  film  probe  was  used  to 
make  RMS  measurements  of  the  turbulent  fluctuations.  The  anemometer  used 
was  a  Thermo  Systems  Inc.  Model  1050.  This  system  had  approximately  a 
200, 000-Hertz  frequency  response  which  depended  upon  the  operating  temperature 
of  the  sensor.  The  RMS  voltage  was  measured  by  a  Hewlett  Packard  true  RMS 
voltmeter.  Autocorrelations  of  the  fluctuating  hot  wire  signal  were  made  with  a 
Honeywell  Correlator  Model  9410.  Probe  positioning  was  accomplished  by  a 
precision  X-Y  table  controlled  by  a  SloSyn  digital  stepping  motor  and  indexer. 

All  data  except  the  autocorrelations  were  recorded  on  a  Hewlett  Packard  Model 
2412B  Data  System. 

4.2. b.  Procedure 

Impact  pressure,  concentration,  and  static  pressure  were  measured  on  the 
jet  centerline.  The  jet  decay  was  found  to  extend  beyond  the  useful  test  section 
of  the  free  jet  wind  tunnel.  At  au  axial  location  of  10  jet  diameters,  radial 
profiles  of  impact  pressure  and  concentration  were  measured.  Turbulence 
RMS  measurements  and  autocorrelations  were  obtained  at  the  same  locations  in 
the  radial  direction  with  the  conical  hot  film  sensor.  It  was  noticed  that  there 
was  a  coupling  between  the  test  cabin  noise  and  the  signal  from  hot  film 
anemometer.  Whenever  the  test  cabin  produced  howling  or  screeching  tones, 
the  hot  film  output  increased  by  a  large  amount.  To  overcome  this  difficulty  a 
constant  arc?,  nozzle  extension  was  fabricated  to  provide  a  shield  for  the  mixing 
region  against  noise  emitted  by  the  test  cabin.  Alter  this  modification  the  radial 

•7 

data  at  —  “  10  were  determined  once  more. 
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4.2. c.  Accuracy 

The  data  signals  were  recorded  at  a  low  level  without  a  signal  conditioning. 

No  amplifiers  or  bridge  balances  were  used  to  keep  errors  introduced  by  these 
devices  out  of  the  measurements.  Pressure  transducers  were  rated  at  0.10% 
full  scale  combined  nonlinearity  and  hysteresis.  The  gas  analyzer  was  found 
to  be  within  1%  accurate  as  determined  by  repeated  calibrations.  The  accuracy 
of  the  Hewlett  Packard  true  RMS  voltmeter  was  found  to  be  within  1.2%  by 
calibration  with  a  precision  signal  generator.  Accuracy  of  the  hot  film 
anemometer  to  indicate  the  random  turbulent  signals  in  conjunction  with  the 
Honeywell  correlator  could  not  be  determined.  However,  steady  state 
calibrations  were  reproducible  within  a  few  percent.  Probe  positioning  was 
calibrated  and  found  to  be  accurate  to  within  .  001  inch  in  6  inches. 

4.3  NORMAL  JET  PENETRATION  AND  MIXING 

The  objective  of  this  experiment  was  to  examine  the  performance  of  an 
injector  designed  to  maximize  turbulence.  The  injector  was  a  Mach  4  contoured 
no  jzle  of  the  same  exit  area  as  the  coaxial  injector.  This  Mach  number  was 
chosen  to  maximize  jet  velocity  at  reasonable  jet  total  pressure.  The  exit  was 
rectangular  with  a  ratio  of  width  to  height  of  10,  with  the  narrow  edge  facing  the 
stream.  This  injector  was  tested  in  a  Mach  2.5  and  3.0  air  stream  to  determine 
the  Mach  number  effects  on  penetration.  Concentration  measurements  were 
taken  in  the  Mach  3.0  stream  only. 

4. 3.  a  Procedure 

Initial  penetration  experiments  were  made  with  the  injector  mounted  on  a 

flat  plate.  Penetration  was  determined  from  Schlieren  pictures  at  Mach  3. 0  and 

2.5.  Then  the  injector  was  mounted  on  the  same  nozzle  extension  that  was  used 

in  the  coaxial  test.  The  flow  field  at  -=■ —  -  10  was  mapped  using  impact 

^eff 

probes  and  gas  analyzer.  Downstream  variation  of  the  peak  helium  concentration 
was  then  determined. 
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Sketch  9.  Schlieren  Photograph  of  a  Normal  Sonic  Jet 
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SECTION  V 

EXPERIMENTAL  RESULTS 

5.1  STRUCTURE  OF  FREE  JET 
5.1.  a.  Velocity  Profiles 

The  Rayleigh  Pitot  formula  was  used  to  determine  the  local  Mach  number 
from  measured  pitot  pressure  and  the  barometric  pressure.  From  the 
measured  total  temperature,  the  static  temperature  and  speed  of  sound  were 
determined  and  then  the  velocity  calculated  from  the  product  of  Mach  number 
times  speed  of  sound.  For  subsonic  flow  the  Mach  number  was  determined 
from  the  isentropic  total  pressure  relation. 

The  absolute  centerline  velocity  profile  for  Mach  numbers  0.8,  1.4,  2.0, 
and  2.5  are  shown  in  Figure  15.  The  similarity  of  the  radial  profiles  is  shown 
in  Figure  16.  These  data  show  no  effect  of  total  temperature  or  Mach  number. 
Also  plotted  is  the  exponential  or  Gaussian  profile  which  seems  to  fit  the  data 
very  well.  Measured  velocity  half  widths  are  shown  in  Figure  17  for  flows  with 
constant  total  enthalpy.  Half  widths  for  flows  with  excess  total  enthalpy  are 
shown  in  Figure  18. 

In  summary  we  can  state  that  Mach  number  and  total  temperature  influenced 
the  mixing  rate  on  the  basis  of  the  centerline  velocity  graph  and  the  half  widths; 
however,  the  similarity  profile  seems  to  be  unaffected. 

5,  l.b.  Total  Temperature 

The  decay  of  the  total  temperature  is  shown  in  Figure  19.  Examination  of 
this  figure  shows  that  the  total  temperature  profiles  do  not  exhibit  a  potential 
core  of  constant  total  temperature  near  the  jet  nozzle.  The  decay  is  continuous 
from  the  combustion  chamber  and  shows  that  turbulence  is  not  relaxed  in  the 
supersonic  nozzle.  The  decay  rate  then  increases  as  the  shear  generated 
turbulence  reaches  the  centerline.  The  similarity  profile  of  excess  total 
enthalpy  Is  shown  in  Figure  20.  These  profiles  taken  downstream,  in  the 
similarity  region  can  also  be  fitted  by  a  Gaussian  profile  but  with  different 
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temperature  half  radius.  The  enthalpy  half  radius  used  to  generate  the 
similarity  profile  is  shown  in  Figure  21.  These  figures  reveal  the  strong 
effect  of  Mach  number  as  increasing  Mach  number  reduced  the  mixing  rate. 

5.  l.c.  Static  Pressure  Defect 


The  most  striking  observation  of  these  experiments  was  the  existence  of  a 
static  pressure  defect  of  considerable  magnitude  for  each  experiment  (Figure  22). 
Again,  we  can  see  the  strong  effects  of  Mach  number.  While  no  absolute 
accuracy  can  be  given  for  these  data  since  no  method  was  available  to  perform 
a  check  on  this  probe,  internal  consistency  from  one  experiment  to  the  next  is 
evident.  These  data  were  reduced  to  turbulence  normal  etross  through 
Equation  151,  Section  IQ.  Assuming  isotropy  on  the  axis,  these  data  were 
converted  to  turbulence  kinetic  energy  K.  These  data  are  shown  in  Figure  23. 
The  turbulence  energy  appears  to  be  expressible  as  a  function  of  the  jet  velocity 
rather  than  the  Mach  number. 


5.  l.d.  Eddy  Viscosity  and  Prandtl  Number 


Tbo  dynamic  viscosity  for  the  compressible  jet  was  determined  from  the 
velocity  decay,  similarity  radial  profile,  and  velocity  half  width  through  Ihe 
formula 


4  JrfZ 


(2035 


which  can  be  derived  from  Equation  188  for  Gaussian  profiles.  Data  were 
obtained  on  the  centerline  from  Mach  numbers  0.2  to  1.?  for  both  the  heated  and 
unheated  jets.  Incompressible  freejet  data  were  obtained  from  Arseyeva 
(Reference  10)  and  treated  the  same  way.  These  results  are  shown  in  Figure  24 
and  are  seen  to  correlate  very  well  with  the  radial  turbulence  pressure.  Shown 
in  the  figure  is  a  line  drawn  through  the  data  having  the  slope  1/ 2.  The  apparent 
trend  of  the  data  deviates  from  this  line  for  the  higher  intensity.  This  trend  is 
attributed  to  compressibility  effects  beginning  with  the  data  points  corresponding 
to  M  =  0.3  and  above. 
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According  t(;  Figure  24  the  incompressible  eddy  viscosity  is  well  represented 
by  the  Prandtl  formula 

€o  "  u*  -^mix  (204) 


where  is  some  fraction  of  the  local  velocity  half  radius  and  that  the 
compressible  viscosity  is  lower  by  an  amount  which  is  Mach  number  dependent. 


Prandtl  numbers  were  calculated  by  the  formula 


^(if)(*'-cP(r«-T*,7r7} 

= _ _ ..VU- 


(205) 


which  was  derived  from  Equation  189  for  Gaussian  profiles.  Large  scatter  was 
found  in  all  the  experiments  with  heated  jets,  probably  as  a  result  of  the  lack  of 
conservation  of  excess  enthalpy  as  noted  in  Figure  H.  The  data  from  the 
constant  total  enthalpy  runs  were  consistent  and  Prandtl  numbers  ranged  from  0.73 
to  0.80. 

o.l.e.  Dissipation 


Dissipation  of  turbulent  kinetic  energy  was  determined  from  the  formula 


-  dK 

(-7T) 


(206) 


which  results  from  Equation  190  for  negligible  diffusion  on  the  centerline. 

These  data  are  shown  in  Figure  25,  which  reveals  decay  law  is  of  the  form  of 
Equation  170.  The  dissipation  length  scale  is  very  nearly  equal  to  the  velocity 
half  radius.  This  scale  is  comparable  with  the  results  of  isotropic  incom¬ 
pressible  turbulence  for  which  the  dissipation  scale  is  of  the  order  of  the 
integral  lateral  scale  of  turbulence.  Considerably  more  scatter  was  obtained 
from  these  data  although  most  experiments  did  exhibit  individually  the  3/  2  power 
dependence,  the  open  squares  ar^  from  experiment  5  and  show  this  dependence 
quite  weli.  No  Mach  number  dependency  was  observed  in  these  data  although 
the  large  scatter  of  the  data  makes  this  result  inconclusive. 
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5. 2  COAXIAL  MIXING 

5. 2.  a  Pitot  Pressures  and  Static  Pressure 

Pitot  pressure  and  static  pressure  radial  profiles  for  one  half  the  mixing 
region  are  shown  in  Figure  26  at  a  downstream  distance  Z/D^  =  10.  The  pitot 
pressure  shows  a  maximum  at  r/D.  =  0.2  and  then  goes  into  a  trough  out  to 
r/D.  =  1.6.  This  trough  is  a  result  of  the  initial  boundary  layer  which  still 
remains.  The  static  pressure  profile  shows  the  existence  of  an  expansion  wave 
in  the  jet  core  which  results  from  the  shook  expansion  system  of  the  conical 
nozzle  used  as  the  fuel  injector.  The  initial  boundary  layer  is  nearly  1. 1  times 
the  jet  exit  diameter  which  is  sufficient  to  drastically  alter  the  effective  boundary 
conditions  for  this  problem. 

5. 2.  b.  Concentration 

The  concentration  profile  shown  in  Figure  27  is  the  average  cf  both  sides 
of  the  mixing  profile.  This  profile  shows  that  the  jet  spreading  has  moved  less 
than  one  half  jet  diameter  radially  at  a  downstream  distance  of  10  jet  diameters. 

5.2. C.  Velocity  Profile 

The  pitot  pressure,  static  pressure,  and  concentration  profiles  were 
reduced  to  Mach  number,  Figure  28.  The  velocity  was  calculated  from  the 
Mach  number  profiles  by  assuming  that  a  constant  total  enthalpy  profile  existed. 
Figure  29.  The.  velocity  profile  shows  the  same  extent  of  spreading  as  the 
concentration  profile.  This  figure,  however,  shows  that  a  large  part  of  the 
velocity  defect  has  been  removed  in  the  boundary  layer. 

5.2. d.  Turbulence 

The  Euierian  integral  time  scale  was  calculated  from  the  autocorrelations. 
The  time  scale  given  by  the  integral  to  infinite  time  delay 

TE  8  /  RrdT  (207) 
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It  is  of  fundamental  importance  in  tbe  present  theory  to  determine  the 
variation  of  the  scale  of  turbulence  across  the  mixing  zone.  By  using  the 
"Taylor  Hypothesis"  that  turbulence  is  locked  into  the  flow  locally;  we  can 
derive  a  longitudinal  scale  of  turbulence  from  the  formula 


A  =  WTe  (208) 

The  results  are  shown  in  Figure  32.  First,  we  see  that  the  scale  is  of  the 
same  order  of  magnitude  as  the  jet  diameter.  Secondly,  v  e  find  that  the  scale 
is  constant  in  the  mixing  layer  and  that  a  separate  scale  ex  sts  in  the  boundary 
layer.  Thus,  this  flow  contains  two  turbulent  shear  layers,  a  new  shear  layer 
due  to  mixing  and  an  old  sheer  layer  from  the  injector  boundary  layer. 


5.3  NORMAL  INJECTION 

5. 3. a.  Optical  Measurement  PenetraUon 

Schlicran  photographs  of  normal  jet  penetration  into  the  air  stream  with 
varying  jet  total  pressure  are  shown  in  Figures  33  and  34  for  free  stream  Mach 
numbers  of  2.5  and  3.0.  As  can  be  seen  by  comparing  these  photographs  with 
the  normal  sonic  jet,  Figure  9,  the  complex  barrel  shock  pattern  does  not 
exist.  Thus,  we  expect  that  more  dissipation  has  occurred  by  turbulence  than 
in  the  normal  sonic  jet.  Penetration  was  taken  to  the  edge  of  the  light  zone 
running  nearly  horizontal  in  the  pictures.  This  definition  corresponds  to  the 
zero  concentration  point  as  determined  by  gas  samples  in  the  next  section. 
Penetration  was  compared  to  the  empirical  correlation  formula  of  Povinelli  and 
Povinelli  (Reference  11),  The  results  are  shown  in  Figure  35.  The  results 
for  the  high  shear  injector  are  correlated  very  well  by  the  equation 
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which  is  Equation  6  from  Reference  6  with  the  slope  increased  by  £5%.  This 
equation  was  obtained  by  Povinelli  using  a  regression  analysis  of  experimental 
data  using  the  dimensionless  groups  determined  by  Vranos  and  Nolan  (Reference 
12).  Since  the  structure  of  this  jet  is  different,  it  is  somewhat  surprising  that 
thia  functional  grouping  works  so  well. 


5. 3.  b.  Concentration  Measurement 


The  concent  ration  profile  normal  to  the  wall  is  shown  in  Figure  36,  10  jet 
diameters  downstream.  This  jet  has  the  same  mass  flow  rate  as  that  in  the 
coaxial  test.  A  very  large  increase  in  mixing  rate  is  observed  *  This  profile 
also  verifies  the  definition  of  penetration  for  the  optical  data.  A  correlation 
of  downstream  decay  of  normal  sonic  jets  was  obtained  by  Henry  and  presented 
in  Reference  7.  Mass  fractions  were  correlated  with  downstream  distance 
divided  by  the  core  length.  Core  length  is  defined  to  be  the  last  downstream 
station  for  which  the  peak  concentration  is  100%.  From  the  Schlieren  photographs 
of  the  jet  structure  and  the  extrapolation  of  concentration  data  back  to  the  100% 
point,  it  is  concluded  that  the  core  length  for  this  jet  is  equal  to  the  half  width  of 
the  nozzle  at  the  exit  with  distance  measured  from  the  jet  centerline, 

_  Width  of  Nozzi# 

X0  = - - -  (210) 


Downstream  decay  of  peak  mass  fraction  was  then  plotted  in  Figure  37  for  a 
comparison  with  the  mixing  rate  of  the  sonic  normal  injector.  The  experimental 
data  are  represented  by  the  formula 


where  the  sonic  jets  were  correlated  in  Reference  5  as 


(211) 


(212) 


These  exponents  show  a  much  faster  mixing  rate,  in  fact  larger,  by  the  square. 
It  also  appears  that  the  mixing  length  is  determined  solely  by  the  size  of  the 
injector  and  not  by  the  dynamics;  however,  more  data  should  be  obtained  to 
verify  this  conclusion.  This  result  would  be  of  great  value  in  engine  design. 
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5.4  SUMMARY  OF  IMPORTANT  EXPERIMENTAL  RESULTS 

This  experimental  work  has  several  important  results  that  apply  directly  to 
both  theoretical  and  practical  design  conclusions. 

a.  Large  defects  in  m»an  static  pressure  are  found  in  compressible 
turbulent  shear  flows.  This  defect  is  interpreted  to  be  a  direct  measurement  of 
the  turbulence  pressure. 

b.  The  experimental  laws  relating  eddy  viscosity  and  dissipation  of 
turbulent  energy  to  the  local  turoulence  energy  and  an  appropriate  scale  can  be 
applied  directly  to  the  formulation  of  an  improved  theoretical  model  of  turbulent 
shear  flow. 

c.  The  scale  of  turbulence  appears  to  be  constant  across  the  shear  layer 
even  though  eddy  viscosity,  dissipation,  and  turbulent  energy  vary  across  the 
shear  layer. 

d.  Turbulent  energy  can  be  frozen  in  a  supersonic  nozzle  and  thus  cause 
performance  losses. 

e.  The  normal  fuel  injector  design  tested  can  apparently  give  short  constant 
mixing  lengths  over  a  wide  range  of  initial  air  velocities,  thus  simplifying 
combustor  design. 
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SECTION  VI 


NUMERICAL  RESULTS  FOR  FREE  JETS 


6. 1  DESCRIPTION  OF  TURBULENCE  PROGRAM  EMPIRICAL  INPUT 

The  equations  for  axisymmetric  mixing  were  coded  for  solution  on  an  IBM 
7094  computer.  The  purpose  of  this  calculation  was  to  show  that  a  coupled  set 
of  equations  could  adequately  describe  the  turbulent  shear  flows.  Many 
empirical  inputs  are  needed  since  the  detailed  turbulent  structure  is  unknown. 
The  following  incompressible  eddy  viscosity  formula  was  used: 

*°  Z  16. 5  (T  K)  (2S 


wmax  “  wmin 


'  dr  h 


(214) 


(215) 


A  constant  mixing  length  seems  appropriate  in  view  of  the  results  for  the 
coaxial  supersonic  mixing  length  scales  (Figure  32). 

The  dynamic  viscosity  was  then  calculated  from  the  equation 

/*T  *  P  €0  {21S1 

The  effects  of  heat  transfer  and  diffusion  on  turbulent  viscosity  have  been 
formulated  by  Qorrsin  in  Reference  >3  but  were  not  applied  in  this  study.  The 
large  Mach  number  effect  was  accounted  for  by  the  following  formulas: 


=  I-  M  (0<  M<  0.5) 


( 2  »’■  i 
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=  i  14  0. 25M}~‘  (0.5  <  M  <  CO  ) 


The  first  formula  is  derived  from  the  definition  of  the  Reynolds  stress  and  the 
assumption  that  the  co variant  density  velocity  fluctuations  are  isentropic  which 
should  be  valid  for  low  Mach  numbers  (Figure  38).  The  second  formula  was 
obtained  by  fitting  the  solution  to  the  data  of  Eggers  (Reference  14),  as  shown 
in  Figure  39.  The  dissipation  function  was  calculated  by  the  formula 


/J.  „n5/2 
3  -  (T  K) 


k  =  2  p 


(219) 


Thus  there  must  be  a  relation  between  the  viscosity  scale  and  the  large  scaie 
eddies,  namely 


'mix 


(220) 


This  value  seems  to  be  physically  plausible.  No  compressibility  or  transport 
effects  were  incorporated  into  the  calculation  of  dissipation. 

The  measured  value  of  Prandtl  number  P  =0,75  was  used  for  the 

r 

calculations  and  the  effective  Lewis  number  was  set  equal  to  unity  for  intuitive 
reasons,  and  then  the  Schmidt  number  became  0. 75  because  of  the  relationship 
between  Schmidt  number  and  Prandtl  and  Lewis  numbers.  The  diffusion  of 
turbulent  kinetic  energy  is  analogous  to  turbulent  heat  and  mass  transfer.  A 
Lewis  number  for  turbulence  defined  below  was  set  equal  to  1, 3  in  the  final 
calculations. 

3  *k 

L»  =  T7f  '22l> 


6.2  SOLUTIONS  FOR  INCOMPRESSIBLE  FLOW 


The  data  of  Laurence  (Reference.  15)  were  used  as  a  comparison  for  the 
low  speed  calculation.  The  lower  curve  in  Figure  39  shows  the  comparison  of 
centerline  velocity  decay.  The  agreement  is  excellent.  This  model  has 
correctly  calculated  the  core  length  as  well  as  the  centerline  decay  in  the 
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transition  and  self-preserving  regions.  Figure  40  shows  a  comparison  of 

measured  and  calculated  turbulence  intensity  in  the  initial  shear  layer  near  the 

z 

end  of  the  potential  core  at  s  3,8.  The  agreement  is  good  but  two 
discrepancies  are  noted.  First,  the  peak  intensity  is  too  low.  This  may  be 
due  to  the  nonisoiropy  of  the  shear  flow,  which  is  not  accounted  for  by  this 
theory.  The  second  discrepancy  is  the  large  difference  in  centerline  turbulence 
level  between  the  calculated  and  measured  values.  The  reason  for  this 
discrepancy  is  not  known.  Perhaps  it  is  due  to  the  assumption  regarding  the 
purely  dissipative  role  of  viscosity  which  may  not  be  valid  at  the  edges  of  the 
turbulent  region.  Some  remarks  in  this  regard  may  be  found  in  Hinze 
(Reference  16).  It  is  interesting  that  the  turbulence  intensity  is  larger  at  the 
outside  edge  of  the  shear  layer  than  the  mean  velocity  itself. 

6.3  SOLUTIONS  FOR  COMPRESSIBLE  FLOWS  WITH  CONSTANT 
TOTAL  TEMPERATURE 

Comparison  of  the  theoretical  calculations  with  the  free  jet  data  for  Mach 
numbers  of  0,8,  1.4,  2.0,  2.5  is  shown  in  Figure  41.  The  agreement  is 
generally  good  except  for  the  Mach  2. 5  case  for  which  the  decay  rate  is  over¬ 
predicted  in  the  transition  region  even  though  the  potential  core  length  is 
predicted  closely.  The  reason  for  this  deviation  is  not  known.  It  is  postulated, 
however,  that  this  deviation  may  be  due  to  the  fact  that  the  turbulence  terms 
were  not  accounted  for  in  the  data  reduction  by  Rayleigh's  pitot  formula.  An 
examination  of  the  Mach  2.0  results  also  shows  a  similar  deviation  but  smaller 
when  the  turbulence  pressure  is  significant, 

6.4  FLOW  WITH  EXCESS  TOTAL  ENTHALPY 

The  low  speed  jets,  M  =  0.8  and  M  =  1.4,  have  heat  transfer  in  the  same 
direction  as  momentum  transfer  while  the  heat  transfer  in  the  high  speed  constant 
total  enthalpy  cases  was  opposite  to  the  momentum  transfer.  Thus  these  cases 
should  present  a  difficult  test  for  the  theoretical  calculation.  The  results  for 
two  cases  are  shown  in  Figures  42  and  43.  Good  agreement  was  obtained.  For 
these  calculations  the  measured  initial  profiles  were  used  as  input.  The 
calculation  was  performed  with  14  input  points  on  the  prof  ile  and  the  calculation 
took  less  than  two  minutes,  a  run  time  characteristic  of  all  these  problems. 

The  calculations  at  M  =2.0  and  M  =  2. 5  did  not  agree  as  well  with  the  data 
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(Figures  44  and  45).  The  data  exhibit  double  structure  with  a  centerline  decay 
that  existed  through  the  nozzle  itself.  These  data  show  that  turbulence  is  not 
relaxed  in  the  nuzzle  so  that  indeed  frozen  turbulence  losses  can  be  expected  in 
supersonic  propulsion  nozzles.  The  present  formulation  of  the  mixing  program 
is  for  only  one  shear  layer,  although  in  principle  the  flow  can  be  divided  into 
different  regions  with  different  scales.  Despite  these  shortcomings  it  is  felt 
that  the  present  formulation  of  the  turbulence  program  is  more  useful  and 
consistent  than  previous  turbulent  calculations  based  on  a  velocity  defect  law  for 
eddy  viscosity.  These  calculations  are  also  useful  for  the  problem  of  calculating 
the  noise  generated  by  jets  since  turbulence  intensity  is  a  necessary  input  to 
Lighthill's  theory  of  noise  from  submerged  jets  (Reference  17).  TLe  final 
Figure,  46,  shows  a  comparison  of  computed  and  measured  relative  turbulence 
intensities  on  the  centerline  for  two  differeut  Mach  numbers.  The  agreement 
is  quite  good  although  the  theory  is  a  little  low  consistent  with  the  low  results 
obtained  at  the  centerline  in  Figure  39  for  the  low  speed  jet  of  Laurence. 
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SECTION  Vll 
CONCLUSIONS 


The  results  of  the  theoretical  and  experimental  study  of  supersonic  mixing 
have  led  to  the  following  conclusions. 

1.  Thrust  losses  due  to  freezing  the  turbulence  kinetic  energy  of  the 
co.nbustor  can  be  appreciable;  however,  these  losses  are  not  as  large  as 
originally  presented  by  Swithenbank  in  Reference  2  and  do  not  constitute  a 
fundamental  performance  limitation  to  the  supersonic  combustion  ramjet  as 
previously  thought. 

2.  Mixing  rate  and  penetration  of  normal  jets  can  be  increased  by 
designing  injectors  to  make  more  energy  available  for  turbulent,  dissipation  as 
was  done  in  *bese  experiments  by  eliminating  the  strong  internal  shock  s.  ucture 
and  maximizing  the  ^t  veloc  .y.  It  is  concluded  that  normal  fuel  injection  offers 
the  designer  advantages  of  the  shortest  and.  most  constant  mixing  lengths  for  a 
very  small  loss  of  potential  engine  performance. 

3.  The  theoretical  model  of  turbulence  based  on  the  turbulent  state  concept 
should  provide  improved  calculations;,  under  widely  different  initial  conditions, 
over  the  previous  theoretical  models. 
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APPENDIX  A 

TABLE  OF  EXPERIMENTAL  CONDITIONS 
FREE  JET  EXPERIMENTS 


E'^periment  No. 

M 

P 

0 

■ 

1 

2.0 

115 

531 

Calibration 

2 

2.0 

120 

960 

Calibration 

3 

2.0 

113 

558 

Bad  Pressure  Data 

4 

2.0 

117 

997 

5 

2.0 

114 

1200 

6 

2.0 

115 

560 

7 

2*  o 

250 

900 

8 

2.5 

253 

1500 

9 

i*  0 

— 

— 

P. ,  Probe  Failed 

10 

if  5 

245 

11 

1.0 

26 

535 

Only  Static  Pressure  Measured 

12 

1.0 

27 

i 

Only  Static  Pressure  Measured 

13 

1.0 

25 

1440 

Only  Static  Pressure  Measured 

14 

1.0 

24 

1100 

15 

1.0 

25 

540 

16 

1.5 

50 

1330 

17 

1.5 

1050 

D  .ta  Acquisition  System 
Malfunctions 

18 

1.5 

50 

550 

19 

1.5 

50 

545 

20 

1*  ^ 

50 

1034 

21 

1.5 

10* 

555 

22 

1.5 

153 

565 
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Experiment  No. 

M 

Po 

a 

Comments 

23 

2.5 

250 

1350 

24 

2.5 

— 

— 

Cooling  Check 

, 

25 

2.5 

— 

— 

Calibrate  PT 

26 

2.5 

— 

— 

Calibrate 

27 

2.0 

95 

1330 

28 

2.0 

50 

1250 

29 

2.0 

120 

560 

PT  Only 

30 

2.9 

120 

560 

PT  Only 

31 

2.0 

115 

1260 

Naval  Ordnance  Laboratory 

Probe  Used 

32 

2.0 

220 

1365 

33 


2.0 


200 


560 


Only  Static  Pressure  Measured 
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Figure  4.  Specific  Impulse  Losses  Due  to  Freezing  Turbulence  Generated 
by  Mixing 
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Figure  5.  Effect  of  Fuel  Injection  Mode  on  Combus  or  Total  Pressure  Losses 


Figure  8.  Free  Jet  Exit  Total  Temperature  Profiles 
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Figure  9.  Schematic  Drawing  of  Instrumentation 
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Figure  12.  Accuracy  of  Exit  Momentum  Flux 
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Figure  17.  Experimental  Velocity  Half  Width;  T  /T 


Stagnation  Tempeature  on  Jet  Centerline 
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Figure  13.  Free  Jet  Centerline  Total  Temperature  Decay 


Experimental  Free  Jet  Centerline 


Viscosity  Correlation  With  Turbulent  Pressure 
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Figure  27.  Profile  of  Helium  Mole  Fraction 
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Integral  length  Scale  of  Turbulence 
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Figure  36.  Profile  of  Helium  Concentration 
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Relative  Turbulence  Intensity 
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Figure  40,  Comparison  of  Computed  and  Measured  Turbulent  Kinetic  Energy 
from  Reference  15 
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